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ABSTRACT 


Muskhelishvili’s  method  of  complex  potentials  has  been  applied  to  the  problems 
of  one  crack,  and  two  diametrically  opposed  (symmetrical)  cracks  emanating 
from  a  circular  hole  of  radius  R,  subjected  to  a  biaxial  load.  The  cracks  of  length 
a,  are  orthogonal  to  the  principal  applied  stress  o  ,  with  transverse  stress 
o  ~  =  Xo  ™ .•  This  work  extends  previous  work  through  the  inclusion  of  linear 
springs  with  spring  constant  k  bridging  the  crack  opening.  Analysis  focussed  on 
the  (normalized)  design  parameters  of  crack  tip  stress  intensity  factor  Fn  and 
crack  mouth  opening  Vn.  Their  dependencies  on  biaxiality  X,  normalized  spring 
stiffness  ka,  and  the  geometry  specified  by  an  =  a/(a  +  R),  were  investigated. 
Interpolation  formulae  with  parameters  depending  on  a„  were  fitted  to  the  high 
and  low  ka  limits  of  F„  and  Vn.  These  provided  a  simple  means  for  calculating 
Fn  and  Vn,  in  most  cases  to  within  a  few  percent  of  the  numerically  calculated 
values.  An  interesting  comparison  of  the  symmetrically  cracked  hole  to  the 
partially  bridged  centre  crack,  showed  that  the  latter  had  a  lower  stress  intensity 
factor  in  all  but  the  very  short  crack  cases. 
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Stress  Intensity  Factors  and  Crack  Mouth 
Openings  for  Bridged  Cracks  Emanating  from 

Circular  Holes 


EXECUTIVE  SUMMARY 


A  frequently  occurring  maintenance  problem  is  the  repair  of  fatigue  cracks 
originating  from  holes  in  structures.  One  method  of  repair  involves  bonding  a 
composite  material  patch  over  the  crack,  however,  the  effectiveness  and 
durability  of  the  repair  must  be  assessed  prior  to  application. 

In  this  work,  Muskhelishvili's  mathematical  method  is  used  to  model  one 
crack,  and  two  diametrically  opposed  symmetrical  cracks,  of  length  a  emanating 
from  a  hole  of  radius  R  in  a  large  thin  plate.  The  patch  is  modelled  by  springs  of 
stiffness  k  acting  between  the  crack  faces  to  oppose  opening.  This  single 
parameter  incorporates  the  moduli  and  thicknesses  of  the  patch,  adhesive  and 
plate.  The  output  parameters  crack  tip  stress  intensity  factor  ( Fn )  and  crack 
mouth  opening  (Vn),  indicate  the  effectiveness  and  durability  respectively  of  a 
proposed  patch  repair.  These  parameters  are  tabulated,  and  simple 
interpolation  formulae  provided  as  functions  of  the  hole  relative  to  crack  size, 
spring  stiffness,  and  biaxiality  of  the  load  applied  to  the  plate. 
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1  INTRODUCTION 


Cracks  in  structures  frequently  originate  from  stress  concentrators  such  as  bolt  holes, 
rather  than  straight  edges  or  from  centres  of  panels.  Stress  analyses  for  these  last 
two  cases  are  presented  by  Tada  et  al. ,  1985  for  a  large  variety  of  loading  geometries 
and  specimen  geometries.  Some  results  for  cracks  emanating  from  circular  holes  are 
also  presented,  but  they  are  not  so  extensive  and  contain  no  information  on  the 
crack  opening  profile.  They  also  do  not  include  springs  bridging  the  crack  faces. 

The  present  work  was  undertaken  to  extend  the  Tada  et  al,  1985  data  by  including 
crack  opening  information.  It  also  treats  the  case  where  crack  opening  is  opposed 
by  linear  springs  between  the  crack  faces  as  shown  in  fig.  1.  It  is  an  analogous 
development  to  the  work  of  Rose,  1987  for  centre  cracks  and  comparison  with  that 
work  leads  to  some  interesting  results. 


The  two  geometries  considered  here  for  cracks  emanating  from  a  circular  hole  in  an 
infinite  plate  are  shown  in  fig.  1  as  (a)  a  single  crack,  and  (b)  two  cracks  emanating 
symmetrically  on  opposite  sides  of  the  hole.  Two  more  cases  were  examined  as 


Figure  1:  The  loading  geometry  and  (a)  a  single  crack  emanating  from  a  hole,  (b) 
the  symmetric  double  crack  case,  (c)  the  edge  crack  and  (d)  the  centre  crack. 
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limiting  behaviours  of  these:  (c)  an  edge  crack  representing  the  large  hole  or  short 
crack  behaviour,  and  (d)  a  centre  crack  in  an  infinite  plate  as  the  small  hole  or  long 
crack  limit. 


Loading  followed  the  Tada  et  al .,  1985  convention  where  uniform  remote  uniaxial 
tension  perpendicular  to  the  crack  is  obtained  for  A  =  0,  whilst  A  =  1  gives  biaxial 
loading.  The  presence  of  the  hole  causes  the  remote  transverse  stress  a^x  to  affect 
the  behaviour,  unlike  the  edge  or  centre  crack  cases  which  are  insensitive  to  this 
stress. 


The  springs  may  serve  to  model  the  effect  of  a  (repair)  patch  applied  over  the  crack, 
leaving  the  hole  clear,  or  fibre-bridging  for  a  hole  in  a  composite  material.  Stress- 
free  boundary  conditions  were  assumed  for  the  hole  surfaces,  so  that  the  hole  may 
not,  for  example,  contain  an  interference-fit  shaft,  rivet  or  other  device  that  applies 
stresses  to  the  hole  surfaces.  The  general  approach  could,  however,  be  extended  to 
these  cases  by  the  introduction  of  appropriate  boundary  conditions. 

The  complex  potential  theory  of  Muskhelishvili,  1953  is  applied  in  the  next  section 
to  obtain  the  appropriate  potentials  for  an  infinite  plate  with  a  circular  hole,  loaded 
in  uniform  uniaxial  or  biaxial  tension.  A  similar  approach  gives  the  potentials  for 
dislocations  inside  or  outside  the  hole.  The  requirement  that  the  crack  faces  be 
stress-free,  allowing  for  the  effects  of  linear  springs  acting  between  the  faces,  leads 
to  an  integral  equation  for  the  distribution  of  dislocations  which  may  be  considered 
to  comprise  the  crack. 

Numerical  solution  of  the  integral  equation  gives  the  dislocation  density,  from  which 
the  crack  opening  profile,  and  the  two  important  parameters  /\tip  and  crack  mouth 
opening  are  obtained. 

Asymptotic  limits  for  Ktl p,  known  analytically,  are  presented  in  section  3,  and  serve 
as  checks  on  the  actual  numerical  results  of  the  next  section.  An  interpolation 
formula  as  a  function  of  spring  stiffness,  depending  parametrically  on  the  loading 
and  crack  geometries,  is  presented  next. 


The  following  discussion  on  crack  mouth  openings  has  fewer  asymptotic  limits  be¬ 
cause  the  crack  mouth  opening  depends  heavily  on  conditions  near  the  hole.  In 
contrast,  Ktip  depends  more  on  stresses  near  the  tip,  which  may  be  remote  from  and 
thus  less  sensitive  to,  the  hole.  The  interpolation  formulae  for  crack  mouth  opening 
suggested  by  these  limits  proved  to  be  unsuitable  for  some  crack  geometries  because 
of  unwanted  divergences.  Alternative  forms  for  the  interpolating  functions  thus  had 
to  be  investigated. 


The  similarity  of  the  symmetrically  cracked  hole  case  to  the  partially  bridged  centre 
crack  of  Rose,  1987  led  to  the  interesting  comparison  in  section  7.  This  indicated 
that  in  general,  except  for  near-tip  bridging,  the  partially  bridged  centre  crack  had 
a  lower  Ktip.  The  geometry  producing  equality  of  both  cases  was  quite  insensitive 
to  the  spring  stiffness. 
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Table  1:  Young’s  modulus  and  other  parameters  for  plane  strain  or  generalized  plane 
stress  conditions. 


parameter 

plane  strain 

plane  stress 

2  fi 

E/(  1  +  u) 

£7(i  +  ") 

E' 

£/(i  -  i/2) 

E 

K 

1 

CO 

(3 -«')/(!  +  *') 

K  +  1 

4(1-1/) 

4/(1  +") 

2  THEORY 


Throughout  the  analytical  part  of  this  work,  the  methods  of  Muskhelishvili,  1953 
have  been  used  to  deduce  the  appropriate  complex  potentials  $  and  together 
with  the  associated  (planar)  stress  and  displacement  fields.  Cracks  were  introduced 
by  expressing  the  crack-opening  profile  in  terms  of  a  dislocation  density  following 
Bilby  and  Eshelby,  1968: 

/* a-\-R 

S(x)  =  2uy(x,  y  — >  0+)  =  /  D(t)dt.  (2.1) 

Jx 

The  crack  mouth  opening  S(R)  will  be  the  maximum  crack  opening  for  the  double 
crack  case  (fig.  1(b)),  but  may  not  be  for  the  single  crack  (fig.  1(b))  which  is  long 
relative  to  the  hole  radius. 

The  procedure  is  to  firstly  calculate  the  stresses  around  a  hole  in  an  infinite  plate 
with  no  cracks,  then  add  the  appropriate  distribution  of  dislocations  while  ensuring 
that  the  hole  surfaces  remain  stress  free,  to  produce  a  crack.  This  imposes  the 
boundary  condition 

&yy(R  <  x  <  R  +  a,y  0+)  — >  0  (2-2) 

and  the  same  for  crxy. 

Prior  to  numerical  solution,  the  resulting  integral  equation  is  supplemented  by  a 
term  incorporating  the  effects  of  linear  springs  which  provide  a  crack  closing  stress 

1  [a+R 

a%(x)  =  E'kuy(x,0+)  =  -E'kj*  D(t)it.  (2.3) 

The  appropriate  modulus  E'  is  given  for  plane  strain  and  generalized  plane  stress 
(Tada  et  al,  1985)  in  table  1. 


2.1  The  Hole 


In  treating  a  plate  with  a  hole  in  it,  a  first  approximation  is  to  ignore  the  hole 
completely,  thus  obtaining  the  (far  field)  potentials  for  a  plate  loaded  in  a  uniform 
remote  stress  state.  These  are 


$“(x,3/)  =  i(C  +  <T~)  +  iC  =  (2.4) 

*“(*.»)  =  -  o  -  <  = 
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Figure  2:  The  regions  S+,S  and  L,  together  with  the  stresses  in  the  rotated 
coordinates. 

with  z  =  x  +  iy  and  C  an  arbitrary  real  constant.  These  may  be  verified  by  substi¬ 
tution  into  the  Muskhelishvili,  1953  relations: 

2fi(ux(x,y)  +  iuv(x,y)y  =  ac^(z)  -  z$(z)  -  ip(z) 

<rxz(x,y)  +  cryy(x,y)  =  2  ($(z)  +  $(z))  =  4Re{<f>(z)}  (2.5) 

(r*x(x,y)- avy(x,y)  +  2iaxy(x,y)  =  -2(z$'(z)  +  T(z)). 

Here  y  denotes  the  shear  modulus,  and  k  depends  on  Poisson’s  ratio  u  as  shown  in 
table  1.  $'(z)  =  The  displacement  of  point  (x,y)  is  ux  parallel  to  the  x  axis, 
and  uy  along  y.  For  the  current  problem,  C  -  0  and  cr“  =  0  in  equation  (2.4). 

The  presence  of  the  hole  modifies  the  potentials  from  (2.4)  to 

Hz)  =  r°(z)  +  $h(z)  (2.6) 

$(z)  =  r(z)  +  $h(z) 

where  $  and  T  are  perturbations  due  to  the  hole.  They  are  analytic  everywhere 
outside  the  hole,  vanishing  at  infinity  faster  than  1/z.  This  rate  of  decay  at  infinity 
is  required  to  ensure  that  there  is  no  net  dislocation  content,  force  or  moment  about 
the  hole. 

The  hole  surface  is  defined  as  L,  points  on  L  denoted  by  t,  and  S+  the  region 

inside  with  S~  outside  the  hole  as  shown  in  fig.  2.  Application  of  the  rotation  of 

coordinates  formulae 

uT  +  iug  =  exp(— i6)(ux  -f  iuy) 

O-rr+Vee  =  Vxx  +  CTyy  (2.7) 

arr  —  &ee  +  2iare  —  exp(— i26)(axx  -  ayy  +  2 icrxy) 

to  equation  (2.5),  together  with  exp(— i'26)  =  z/z,  allows  the  stress-free  hole  surface 
boundary  condition  to  be  written 

CTrr(t)  +  iarg(t)  =  0  =  <f>(z)  +  $(I)  -  z$'(z)  -  {z / z)W{ z)  (2.8) 

as  z  — >  t  from  S~ . 

The  idea  now,  is  that  we  seek  to  re-express  this  in  the  form 

0  =  $(z  — >  t~)  +  H(z*  — ►  t+)  (2.9) 
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where  the  function  f i(z*)  need  only  be  defined  for  z*  £  S+  being  an  as  yet  unspecified 
function  of  ^  G  5".  The  only  requirement  on  2*  is  that,  as  z  approaches  t  £  L  from 
S+,  then  approaches  t  from  S~.  fl(;z*)  will  be  expressed  in  terms  of  $(z)  and 
\J/(z)  such  that  (2.8)  is  replaced  by  (2.9)  in  the  z  — ■»  t  limit. 


If  this  can  be  done,  then  the  function  F(z)  defined  by 


F(z) 


$(z)  ,2  G  S~ 
— f l(z),z  6  S+ 


(2.10) 


is  analytic  across  the  boundary  L.  Its  behaviour  for  z  6  S~  is  dictated  by  that 
for  $~(z).  The  behaviour  at  00,  and  any  poles  in  S~  are  thus  specified.  For  the 
behaviour  in  S+ ,  we  look  at  the  singularities  determined  by  the  form  of  Cl(z*): 
these  may  occur  at  z*  — >  0,  and  at  points  determined  once  more,  by  the  known 
singularities  in  $(z)  and  $(2). 


F{z)  is  therefore  a  function  analytic  in  the  whole  plane  except  at  known  poles  with 
known  coefficients,  and  a  prescribed  behaviour  at  00.  It  must  be  of  the  form 


no -EE  77 

j  n=  1 


+  ^2  dkzk . 


(2.11) 


This  gives  §(z)  directly  for  z  £  S  ,  and  ^(^)  may  be  found  by  inverting  the 
expression  for  Q(z*). 


In  the  present  problem,  a  suitable  expression  is 

fi(z-)  =  $(R2/z*)  -  (R2/z*)W(R2/z*)  -  (R2/z*2P(R2/z*)  (2.12) 

with  2*  =  R2/z  as  suggested  by  List,  1969.  Equation  (2.8)  then  becomes 

o' 'tt(z)  +  ivre(z)  =  $(2)  +  tt(R2 /z)  +  ~  -J  ®(«).  (2.13) 

As  z  — >  t  £  L  from  S~,  (R2 /z)  — >  t  from  S+  and  the  last  term  above  vanishes.  The 
left  side  becomes  zero  due  to  the  stress-free  boundary  condition  (2.8).  This  equation 
reduces,  on  L  as  required,  to  (2.9).  The  analytic  continuation  arguments  above  lead 
to  the  function  F(z)  given  by  (2.11),  with  singularities  prescribed  by  equations 
(2.4),  (2.6),  and  (2.12).  From  (2.6)  and  (2.4),  N  =  0  and  d0  =  \(a ~  +  <r~).  The 
other  singularity  in  F(z)  arises  as  z  — *•  0  in  the  last  term  of  (2.12):  f l(z  — >  0)  — + 
—  (R2  /  z2)'$>{\z\  — >  00).  Finally, 

f w = -  o + j(<c + o  (2-14) 


Equation  (2.12)  may  be  rewritten  to  make  ^(.z)  the  subject: 

V(z)  =  ( R2/z2)$(z )  -  (R2/z)<F{z)  -  ( R2/z2)T}(R2/z ). 


(2.15) 


Making  use  of  (2.11)  and  (2.14),  and  noting  that  if  z  6  S  then  R2/z  €  S+,  we 
obtain  both  potentials  as 


*-(2) 


(2.16) 


*-(*) 


1-1  R2  3/<4  \  „  ( 1  R 2  3fl4\ 

V  2  +  2z2  2 z4  J  +  °yv  \2  +  2z2  +  2z4  )  ' 
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2.2  Dislocation  Outside  the  Hole 


In  addition  to  the  potentials  developed  in  the  previous  section,  we  need  the  potentials 
for  a  dislocation  outside  the  hole  as  in  fig.  3,  and  also  those  for  one  inside  the  hole, 
treated  in  the  next  section.  We  begin  with  the  potentials  for  a  dislocation  at  the 
origin  in  an  infinite  plate: 


®°(*) 


V°(z) 


—ft  iBx  —  j By  —i  A 

7r(  1  +  K)  Z  Z 

ft  iBx  +  By  _  i A 

7T  (1  +  k)  Z  Z 


(2.17) 


where  A  =  and  B  =  Bx  +  iBy  denotes  the  Burgers’  vector.  List,  1969 

appears  to  have  a  sign  wrong  in  his  equation  2.7  which  should  have  —  C  —  D.  The 
above  are  correct  by  substitution  into  equations  (2.5).  Application  of  the  translation 
of  coordinates  formulae 


ZB  =  ZA  ZBA 

$b(zb)  =  $a{za) 

Tb(2b)  =  ^a{za)  +^ba$'b(zb) 

gives  the  potentials  due  to  a  dislocation  at  b  as 


*D(z) 


—i  A 

z  —  b 

i  A  iXb 

z  —  b  (z  —  b)2 


(2.18) 


(2.19) 


The  hole  modifies  these  due  to  the  stress-free  hole  surface  in  much  the  same  way 
as  for  the  remote  applied  stress  case.  Equations  (2.6)  to  (2.11)  still  apply  with  <f>°° 
and  T00  replaced  by  $D  and  TD  respectively.  The  singularities  in  F(z)  as  given  by 
(2.11)  will  differ,  so  from  there  a  separate  treatment  is  required. 
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Singularities  in  $D  and  ’f0  produce 

Qsingular/  X  =  ^  _  R*  ~iX  ,  &  ~^b 

(R2/z)  —  b  *  ((B?/z)-b)2  z2(R2/z)-b  z2  ((B?/z)-iy 

(2.20) 

Using  d  =  R2/b ,  expanding  using  partial  fractions,  and  collecting  pole-type  terms, 
F(z )  becomes 


.  —i X  i X  iX  iXd(b-d) 

F^  =  —i  +  —J-7  +  T(hw 


(2.21) 


Equation  (2.15)  becomes,  using  (2.21)  in  (2.11), 

—2iXR2  iXR2(2z-d )  iXR2(2z-b) 

^w)  —  „3  ~2  (~  J\2  ~2(  ~  U\2 


\z  -  dy 


\z  -  by 


(2.22) 


iXd2(b  —  d)(3z  —  d)  iX  iXd  iXb  iX{b  —  d) 

d~  n  /  i\o  U  /  T\  /  7  \  /  TV1 


l(*  -  dy 


z(z  —  d)  z(z  —  b)  (z  —  b)2 


This  may  be  simplified  by  partial  fractions,  producing 

iX  iXb  iX  iX  2iXR2 
=  ~  (z  -  b)2  "  z  -  d  +  7  7~ 

+i_(A(6  -d)-  X(b  -  d))  +  -  1(6  -  d)) 

2iXd(b  —  d) 

+  (z-df  ■ 


(2.23) 


2.3  Dislocation  Inside  the  Hole 

The  potentials  are  now  sought  for  a  dislocation  located  inside  the  hole,  for  conve¬ 
nience  at  the  origin.  These  are  then  combined  with  those  from  the  previous  section 
to  render  no  net  dislocation  content  far  from  the  hole.  The  procedure  is  the  same  as 
before,  but  and  D  are  replaced  by  $°  and  'k0  from  equation  (2.17).  The  only 
singularity  in  F{z )  by  (2.11)  is  as  z  — >  0  in  the  O(z)  term:  S7(z  — »  0)  — >  y.  This  is 
from  the  —{R2/z2)^{R2jz)  term.  Equation  (2.21)  is  replaced  by 


F\z)  = 

and  substitution  into  (2.15)  produces 


*°(*), 


(2.24) 


z  z 6 


(2.25) 


This  has  been  superscripted  “Oh”  to  distinguish  it  from  \P°  in  equation  (2.17). 


2.4  Stresses  Along  y  =  0 


Knowing  the  potentials  for  all  three  cases,  the  stresses  ayy(x,y  =  0)  they  produce 
need  to  be  calculated.  These  stresses  enter  into  the  crack- defining  condition  (2.2) 
and  produce  the  integral  equation  for  the  dislocation  density  and  hence  crack  profile. 


We  put  z  —  x  in  (2.16),  (2.21),  and  (2.23)  to  (2.25),  and  substitute  these  into  (2.5) 
to  obtain  cryy(x,  0).  Furthermore,  b  is  taken  to  be  real  so  that  d  —  R2/b  is  also 
real,  and  A  =  iBn ,  that  is  Burgers’  vector  is  purely  imaginary  ( iBy  only).  These 
simplifications  make  3>(x)  and  \k(a:)  real,  where  the  y  =  0  has  been  omitted  for 
brevity.  Equations  (2.5)  indicate  that  crxy(x )  =  0  as  required,  and 

cryy(x)  —  2<h(:r)  +  z$'(x)  +  'I'(x).  (2.26) 

Using  this  equation  and  the  above  simplifications,  the  following  results  are  obtained 
where  Bn  = 


1.  Uniform  remote  (biaxial)  tension: 

o„(x>R)  =  (l  +  \{RI*f  +  \w*y)  (2.27) 

+C  Q(*/*)2  -  |{JR/x)-)  . 

2.  Dislocation  located  at  (real)  b  >  R: 

.  (  1  1  1  b  —  d  R2 

(Tyy(x  >  R)  =  2 Bn< - -H - -H - -  I — -  (2.28) 

| x  —  b  x  x  —  d  x *  x J 

(t  —  d )2  (r  —  d )3 

3.  Dislocation  with  Burgers’  vector  B'  =  iB'y  with  B'n  =  located  at  b  =  0 

inside  the  hole: 

ayy{x  >R)  =  2B'n  +  (2  ‘29) 

Burgers’  vector  is  dashed  here  to  distinguish  it  from  the  previous  item:  later 
we  will  set  B'n  =  —Bn. 

4.  Putting  R  —*  0  implies  d  — >  0  and  only  the  first  term  remains  in  equation 
(2.28).  This  is  the  stress  due  to  a  dislocation  at  b  in  an  infinite  plate. 

O  D 

^(*)  -  ^rb  (2.30) 

which  could  equally  well  have  been  derived  directly  from  equations  (2.26)  and 
(2.19). 

5.  The  other  limiting  case  occurs  as  R  — »•  oo,  when  the  situation  becomes  a 
dislocation  located  at  b  from  an  edge  at  x  —  0.  This  is  shown  in  fig.  4,  and 
treated  in  the  next  section. 


2.5  Dislocation  Near  an  Edge 

The  potentials  for  a  dislocation  near  an  edge,  as  shown  in  fig.  4,  may  be  derived 
in  a  similar  procedure  to  that  for  the  hole  case.  The  free-boundary  condition  still 


Figure  4:  A  dislocation  of  Burgers’  vector  B  —  Bx  +  iBy  located  at  point  b  in  a  body 
occupying  region  S~  with  its  edge  L  along  x  —  0.  Along  the  edge,  the  boundary 
condition  is  cxx  +  icrxy  —  0.  An  image  dislocation  is  located  in  A"*"  at  —b. 

applies,  but  now  to  the  edge  rather  than  the  hole  surface  of  the  previous  case.  From 
equations  (2.5),  with  the  edge  along  x  =  0, 

Vxx(x,y)  +  i(?xy(x,y)  ~  $(z)  +  $(z)  —  z$'(z) —  y>(z)  (2.31) 

->0  asi->0+. 

For  the  edge  case,  equation  (2.12)  is  replaced  by 

ft(z  e  S+)  =  $(-z)  -  zW(-z)  -  V(-z)  (2.32) 

so  that  _ 

axx(x,  y)  +  icrxy(x ,  y)  =  $(z)  +  -  (z  +  z)$'(z).  (2.33) 


As  z  ^  t  £  L,  axx+iaxy  must  vanish,  and  z  -z  &s  z  becomes  purely  imaginary,  so 
that  once  again,  $_(t)+f2+(t)  =  0.  Using  the  same  analytic  continuation  arguments 
as  before,  and  noting  the  prescribed  singularities  from  equations  (2.19)  and  (2.32), 
we  obtain  F(z)  as 

Fu)  =  zib.  +  AL  +  2A+3-.  (2.34) 

[  ’  Z-h  z  +  b  (2  +  f>)2 

Equation  (2.32)  may  be  rewritten  as 

tf(z)  =  $(*)  +  Z&{z)  -  Ti(-z),  (2.35) 

from  which 

t  /  N  iX  iXb  iX  t  iXb  iX(b  +  b)  (  2iXb{b  +  b)  (f) 

^^)  =  ^_(z_fe)2_^|+(2  +  6)2"  (2  +  6)2  +  (2  +  £)3 


These  potentials,  with  the  same  simplifying  assumptions  as  before,  produce  ayy 
stresses,  on  the  positive  real  axis,  of 


/_ L _ 1_  _  2b  4 b2  \ 

\x  -  b  x  +  b  (x  +  b)2  ^  (x  +  b)3  J 


(2.37) 


=  2  B 


n 


8b2  x 

(x  —  b)(x  +  b)3 
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Figure  5:  Crack  elements  used  for  the  cases  of  fig.  1.  In  each  case  except  the  edge 
(c),  the  element  actually  consists  of  the  superposition  of  two  dislocations  and  the 
stresses  they  produce. 

2.6  Integral  Equation  for  the  Dislocation  Density 


The  crack  opening  profile  is  now  calculated  by  representing  the  crack  as  a  distribu¬ 
tion  of  dislocations  along  the  x  axis  according  to  equation  (2.1).  The  stress  cryv(x ) 
is  then  equal  to  the  sum  of  all  the  stresses  caused  by  these  dislocations,  that  due 
to  the  (remote)  applied  load,  and  the  term  (2.3)  due  to  the  bridging  springs.  This 
total  must  be  zero  for  each  x  within  the  crack,  as  the  crack  surfaces  are  stress-free 
according  to  (2.2). 

The  form  of  the  dislocation  density  depends  on  the  geometry  of  the  problem  being 
solved.  For  the  cases  (except  the  edge  crack)  being  examined  here,  it  is  simpler  to 
take  the  crack  elements  as  dislocation  pairs  as  shown  in  fig.  5  rather  than  individual 
dislocations.  This  permits  automatic  satisfaction  of  the  no  net  dislocation  content 
and  symmetry  constraints. 


Taking  the  single  crack  from  a  hole  case  (a)  in  fig.  5,  the  crack  element  causes  a  stress 
for  R  <  x  <  R .+  a  given  by  equation  (2.28)  minus  (2.29)  with  Bn  =  D(b)db  = 
B'n.  The  total  stress  at  x  due  to  the  crack  is  then  the  integral  of  this  with  respect 
to  6: 


(x) 


2 


rR+a  /  1 

Jr  \x  — 


tt(1  +  k)  Jr  \x  —  b 
d(d  —  b  +  d(l  —  d/b)j 
( x  —  d )3 


D{b)db 


b  —  d  d  —  b 

'  o  f 


(.r  —  dy 


(2.38) 
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where  d  —  R2/b. 


The  integral  equation  for  the  dislocation  density  then  results  by  combining  this,  the 
term  due  to  springs  (2.3)  and  that  due  to  the  applied  load  (2.27)  in  (2.2): 

</(*)=  0  =  — f(x,  6;  R)D(b)db  -  l-E'k  j^‘  D(b)ib  (2.39) 

+<,  i1  +  \WX)2  +  |(^/^p)  +'<C  (\w*)2  -  \w*)')  ■ 

For  the  single  crack  from  a  hole,  f(x,  b]  R )  is 

A  1  1  b  —  d  d-b  d(d-b  +  d{l  -  d/b)) 

F  b ’ R)  ~  ^ZTh  ~  ^Z7d  +  x2  +  (x  _  d )2  ~  {x-df  (  '  0) 

whilst  for  the  symmetric  double  crack  it  is  obtained  from  (2.28)  and  the  same  sub¬ 
tracted  after  replacing  b  and  d  by  —  b  and  —  d  respectively: 

_  1  1  1.1.2 (b-d) 


+a7y  ( 1  + 


Fhl(x,  b ;  R)  = - -  — - -  H - ~2  b 

x  —  b  x  -  d  xz  (x  -  d)z 


(2.40) 


.  _  .  1  1  1  1  2  (b-d) 

Fh2(x,b-,R)  =  - - - — r - + - 2 

x  —  o  x  +  b  x  —  d  x  +  d  xz 

+(d  ~  ^  (( x-d )2  +  (. x  +  dy ) 


(2.41) 


-d{d-b  +  d(l-d/b)) 


V  ‘  ^  ,  J'\(x-d) s  (x  +  dfj' 

The  centre  and  edge  crack  cases  are  obtained  by  setting  R  =  0  in  (2.39),  reducing 
the  remote  loading  term  to  cr~.  The  integration  ranges  become  0  to  a  and  x  to  a. 
f(x,  b ;  R)  becomes  a  function  of  (x,  b )  only,  obtained  from  (2.37)  as 


for  an  edge  crack,  and 


£,  ln  8b2x 

^  ( x-b)(x  +  b )3 

(2.42) 

f^=x1-b-x\b 

(2.43) 

from  (2.30)  for  the  centre  crack  case. 

2.7  Normalization 


Prior  to  discretization  and  numerical  solution,  equation  (2.39)  must  be  normalized, 
with  a  natural  choice  for  the  stress  and  length  (fig.  1)  scales  being 

O  _  _F  _ 1  7  _  D  fr,  A  A\ 


Sn  —  .  —  .  and  ln  —  R. 

7r(l  +  k)  47T 


(2.44) 


For  edge  and  centre  cracks,  ln  =  a.  The  following  normalized  variables  are  then 
defined: 

Sx=a™/Sn  S™h(X)=a™h(x)/Sn 
X=x/ln  B=b/ln  D=d/ln  =  l/B.  1  ' 

The  stress  function  is  normalized  as F(X,B)  =  lnf(x,b]R )  while  the  dislocation 
density  is  already  normalized  because,  from  (2.1),  D(b)  =  —  =  -d^^.  The 

length  of  the  crack  relative  to  the  hole  is  specified  by 


an  =  a/(R  +  a) 


(2.46) 
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Table  2:  The  normalized  limits,  Bi  and  Br  for  the  first  integral  in  equation  (2.39) 
as  required  by  the  four  crack  geometries  of  fig.  1. 


Case 

B, 

Br 

In 

hole  (a),(b) 
edge  (c),  centre  (d) 

1 

0 

K+a.  _  1 

R  —  l-a„ 

1 

R 

a 

so  that  short  cracks  imply  an  — >  0  whilst  long  cracks  have  an  — >  1. 

The  limits  on  the  first  integral  in  (2.39)  become  Bi  and  Br,  given  in  table  2  for  the 
four  cases. 

After  normalization,  equation  (2.39)  becomes 

fBr  F(X,  B)D(B)dB  -  2,ir kln  D(B)dB  =  S™h(X)  (2.47) 

J  B[  J  X 

where,  for  the  hole  cases 

Syy  C X )  =  Sy(l  +  —  +  2X4  )  +  S*  (2X2  ~  2X1)  '  (2-48) 

For  the  edge  and  centre  cracks,  S™(X)  —  Sy  only.  The  main  results  required  from 
(2.47)  are  the  crack  tip  stress  intensity  factor  Ktip  and  crack  mouth  opening  ln8n{B{). 
These  are 


tftip  =  —  \\m_{y]2i7lnyj Br  -  BD(B)} 
=  lim  {y/2irlnyjB  -  BrSnSyy(B)} 

B — yBT 

6n(Bi)  =  fBr  D(B)dB. 

J  Bi 


(2.49) 

(2.50) 


The  first  of  these  demonstrates  the  singularity  expected  for  D(B )  for  any 

crack  with  a  non-zero  Kt ;p.  To  overcome  this  problem,  and  increase  the  accuracy  of 
the  numerical  work  by  decreasing  the  discretization  intervals  as  B  — >  Br  where  D(B ) 
changes  most  rapidly,  a  variable  transformation  is  made.  Following  the  procedure 
of  Rose,  1987, 

B  =  Bt  sin (t)  dB  =  Bt  cos(t)df 
v  _  d  n/o\  _  Q(t) 


X  =  Br  sin(u)  D[B)  =  ,  (2.51) 

Br  COS (t) 

ti  =  sin  ~l(Bi/ Br)  tr  =  7t/2 
and  the  integral  equation  becomes 

F(BT  sin(u),  Br  sm(tf)Q(t)dt  +  2xkln  J  Q(t)dt  =  S™h(Br  sin(u)).  (2.52) 

The  interval  [ti,tT]  is  split  into  N  uniform  intervals  of  width  St  with  midpoints 
t(j)  —  U  +  ( j  —  \)8>t  as  in  fig.  6.  The  first  integral  for  u  =  t(j)  can  then  be  split  up 
as  I(j)  =  Eili  SI(jti)  where 

Q(t)F{t(j),t)dt.  (2,53) 
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12  3 

•  X  !  X  1  X 
sin'7  (bjb) 
X=Bl 


j  N-2  N-l  N 

- 1 . x-b - 1  ><  1  X  1  X  • 

t(j-  V  t(j)  \(j+ 1)  tr=B 

X(j)=Brsm(t(j))  X=B 


Figure  6:  The  discretization  of  the  interval  [f/,  tT]  showing  N  intervals,  the  jth  having 
boundaries  t(j  —  |)  and  t(j  +  ^). 


Two  approximations  are  now  made:  firstly  Q(t)  is  taken  as  a  constant,  Q(i)  — 
Q(t(i)),  over  the  interval  of  integration,  and  secondly  the  trapezoidal  rule  is  used  to 
evaluate  the  remaining  integral.  These  give 

=  Q(i)6j{F(t(j)Ai  -  i))  +  +  |))}-  (2.54) 


Particular  care  is  required  when  j  =  i  because,  in  this  case,  the  singularity  in 
F(t(j),t )  needs  to  be  integrated  over.  It  arises  from  the  term  in  F(X,B).  The 
singular  part  is  thus 


si‘,n^“UJ)  =  QU)  l 


r*  0)+f 

1 

kj)-%  BT\ 

sin  (t(j))  -  sin(f)} 

d  t. 


(2.55) 


Changing  the  variable  t  =  t(j)  +  ^Stv  and  expanding  sin(t)  in  a  Taylor’s  series  about 


sin(t)  ~ 


we  obtain 


sin  ( t(j ))  +  (t  -  t(j ))  cos  ( t(j ))  -  |(#  -  t(j))2  sin  (t(j)) 

~  til))*  cos  +---i 


S2t  2 

— v2  + 
24 


-l 


(2.56) 


(2.57) 


The  integrand  may  be  expanded  using  the  series  for  1/(1  —8)  and  the  integral  of 
-  vanishes  by  the  principal  value.  The  integral  of  the  (third)  term  proportional  to 
v  vanishes  too,  by  symmetry.  Just  the  second  term  is  left,  together  with  terms  of 
order  that  is 


<f/smgular(i,j) 


sin 

cos 


+  0(8?). 


(2.58) 


The  same  result  would  have  been  obtained  directly  from  (2.54),  that  is  the  singularity 
is  taken  care  of  by  application  of  that  formula  for  j  =  i. 


Turning  to  the  next  integral  in  (2.52),  the  term  due  to  springs  is 

7sp(i)  =  2tt kln  fr  Q(t)dt  =  2tt kln6n(j).  (2.59) 

Jt(j) 

The  trapezoidal  rule  is  again  used,  but  extra  care  as  t  — >  tr  gives  expressions  for 
Q(t),D(t)  and  8n(t )  that  extrapolate  to  t  =  tT ,  the  one  for  Q(t)  being  useful  later 
to  calculate  Afiip. 
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If,  in  addition  to  (2.51)  the  change  e  =  (tt/2 )  —  t  is  made,  then  the  B  — >  Br  expansion 
of  6n(B)  becomes 


fin{B —>  Br)  —  (a0  +  ot\ (Bt  —  B )  +  •  •  •  ^  \j Bt  —  B 

=  \J BT/2  e (olq  +  ( Btq.i /2  —  ao/24)e2  +  ■  •  • ) 


(2.60) 


Q(e)  =  ~  =  yt/2(Qo(l  -  e2/ 8)  +  3alJBre2/2) .  (2.61) 

Substituting  in  e  =  6tf 2  and  3^</2  gives  equations  for  Q(N )  and  $(./V  —  1)  respec¬ 
tively,  which  may  be  solved  for  <a0  and  : 


a0  = 


ai  = 


[Y  f9Q(Ar)  -  Q(N  —  1) 

Br  8 

-1  [Y  \Q{N)-Q(N  -1) 


(2.62) 


WrMBr[  6 2  J- 

These  may  be  substituted  back  into  (2.60),  leading  to 

K(N)  =  ^(13Q(N)-Q(N~1)) 

S*(N-1)  =  |j(27£?(AT)  4-  9Q(N-lj). 
The  trapezoidal  rule,  in  the  form 

f>n{j)  =  6n(j  +  1)  +  lz{QU)  +  QO  +  1)) 


(2.63) 


(2.64) 


is  then  used  in  (2.59)  to  express  Isp(j)  as  a  matrix  multiplying  the  array  Q(i): 


Pp(j)  =  2n 


(2.65) 


Msp(i,0  = 


12 

24 

•••  24 

21 

27 

0 

12 

•••  24 

21 

27 

0 

0 

•••  12 

21 

27 

0 

0 

...  o 

9 

27 

0 

0 

...  o 

-1 

13 

Setting 


cu)  =  s~h(xu)), 


a  matrix  equation  is  obtained  for  Q(i): 

ZM(j,i).Q(i)  =  C(j) 


(2.66) 


(2.67) 


—  2irkln~Msp(j,i). 
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2.8  Normalized  Ktip  and  Crack  Mouth  Opening 


Returning  to  (2.49)  and  substituting  in  the  derivative  of  (2.60)  for  D(B),  the  crack 
tip  stress  intensity  factor  KtiP  may  be  given  in  a  similar  normalized  form  to  that  of 
Rose,  1987,  that  is 


vTv\f™Fn{ka\  an\  A) 


(2.68) 


Fn(ka,  ctn,  A) 


I  in  /  *  \ 

2a  UJ“0 


with  the  last  equality  for  the  hole  cases  only.  For  edge  and  centre  cracks,  the 
factor  is  replaced  by  l/y/2.  Equation  (2.62)  provided  an  expression  for  a0. 


Both  the  normalization  for  Ktiv  and  the  dimensionless  spring  constant  ka,  could  have 
been  selected  differently,  for  example  the  a  in  a^-JFa  could  justifiably  be  replaced 
by  2 R  +  a  for  the  single  crack  from  a  hole,  and  R  +  a  for  the  symmetric  double 
crack.  The  spring  constant  could  have  been  JcR  but  the  chosen  normalizations  lead 
most  readily  to  the  asymptotic  limits  of  small  hole,  large  hole,  weak  springs  and  stiff 
springs. 


The  crack  mouth  opening  S(b[)  may  also  be  normalized  as 

(4cr°°a\ 

-^-jK(k;cn;A)  (2.69) 

V4ka]a,tX)  =  k  (£  s4Bl)  =  (i) 

The  numerical  integration,  (2.64),  does  not  give  6n(Bi)  because  Sn(j  =  1)  corre¬ 
sponds  to  6n(ti  +  St/ 2)  from  fig.  (6).  We  choose,  in  preference  to  the  trapezoidal 
rule, 

*„(*,)  =  Ml) +  |q(1)  (2.70) 

to  be  substituted  into  the  previous  equation.  Again,  edge  and  centre  cracks  do  not 
have  the  factor,  and  the  dependence  on  A  is  trivial  for  it  has  no  effect  in  these 
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3  ASYMPTOTIC  LIMITS:  Ktip 


As  indicated  previously,  there  are  several  asymptotic  limiting  cases  for  the  cracks 
emanating  from  a  hole,  both  in  terms  of  the  geometry,  R  relative  to  a,  specified  by 
an,  and  the  spring  stiffness  indicated  by  ka.  In  the  following  sections,  stiff  springs 
will  be  treated  first,  large  holes  or  short  cracks  next,  small  holes  or  long  cracks,  and 
finally  weak  springs. 

3.1  Stiff  Springs 


Here,  the  springs  are  so  stiff  that  k  1  is  small  relative  to  both  a  and  R ,  so  that  from 
Rose,  1987,  Kt\P  — >  ayy/Vk-  From  equations  (2.27)  (and  2.45), 


so  that 

Fn(ka  -4  oo;  an;  A)  ->  -^==  ^1  +  ^(1  +  A)(l  -  anf  +  ^(1  -  A)(l  -  an)4j  .  (3.2) 

3.2  Large  R  or  Small  a 


The  simplest  limit  is  R  — >  oo  relative  to  a,  or  an  — ►  0.  In  this  case,  both  hole 
configurations  approximate  to  an  edge  crack  of  length  a,  loaded  by  a  stress  given  by 
the  R/x  — +  1  limit  of  (2.27): 


1 

TT 

si 

s> 

b 

'  3<^“  -  <  =  (3  -  A )SnS„. 

(3.3) 

In  this  case, 

F n^kci,  dn 

-.0;A)-*(3-A  )F^(ka). 

(3.4) 

For  stiff  springs,  from  (3.2),  the  limit  is 

Fn(ka  — >  oo;  an  — » 

0;  A)  — >-i=(3-A-(7-5AK). 

(3.5) 

3.3  Small  R  or  large  a 


As  the  crack  becomes  much  longer  than  the  hole  radius,  an  — >  1,  both  hole  cases 
approach  (different)  centre  crack  geometries,  but  the  asymptotes  may  be  taken  in 
different  ways  for  each  case.  Although  they  limit  to  the  same  in  each  case,  they 
produce  different  forms  for  (2.68)  and  (2.69). 

One  procedure  is  to  take  the  correspondences  shown  in  fig.  7  (a)  and  (b)  for  the 
single  and  double  cracks  respectively,  simple  in  terms  of  variable  substitutions,  but 
not  as  accurate  as  cases  (c)  and  (d),  which  should  hold  further  from  the  an  — >  1 
limit. 

Taking  fig.  7  (c),  the  limit  for  Ktip  becomes  the  fully  reinforced  ( c/a  — >  0)  treatment 
of  Rose,  1987,  where  effects  due  to  the  hole  at  the  other  end  of  the  crack,  including 
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Figure  7:  Possible  correspondences  between  the  hole  configurations  and  an  — >  1  long 
crack  limits.  In  cases  (a)  and  (b),  a  in  the  centre  crack  case  becomes  a/2  and  a  for 
the  single  and  double  crack  cases  respectively.  In  (c),  a  is  replaced  by  (2 R  +  a)/2, 
and  for  (d),  a  becomes  R  +  a. 

the  stress  concentration  of  (2.27)  and  the  lack  of  springs  across  the  hole,  become 
negligible. 

For  no  springs  and  correspondence  (c),  Kt ip  =  o-™\Jtt(2R  +  a)/ 2  whereupon 

c„( 0;  an  -  1;  \)  -*  1  «  /|(l  +  (1  -  aj).  (3.6) 

As  the  springs  become  stiff,  we  take  an  — »  1  in  (3.2): 

Fn(ka  ->  oo;a„  ->  1;  A)  -  -^(l  +  i(l  +  A)(l  -  an)2)  (3.7) 

which  is  the  same  as  Rose,  1987  for  A  =  0. 

The  last  equation  would  still  hold  for  the  double  crack,  case  (d)  in  fig.  7,  but 
equation  (3.6)  wouldn’t  because  Kt iP  is  now  given  by  (x™yJx(R  +  a).  Instead, 

Fn{Q-,an  1;  A)  -►  w  1  +  ^(1  -  an).  (3.8) 

V  ®7l  2 


3.4  Weak  Springs,  ka  — >  0  Limit 

Single  and  double  crack  hole  cases  are  presented  by  Tada  et  a/.,  1985  for  no  springs, 
ka  =  0,  where  s  is  here  represented  by  an,  and  a  by  <r^.  The  other  variables,  R,  a 
and  A  are  the  same. 

The  single  crack  case,  from  page  19.2  in  Tada  et  al. ,  1985,  is  presented  as 
Fn(0;an;A)  =  (1  —  A)Fo(a„)  +  AFi(an) 
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(3.9) 


F0(an)  =  (l  +  0.2(1  -an)  +  0.3(1  -a^F^) 

Fj(an)  =  2.243  -  2.640an  +  1.352a*  -  0.248a®. 

In  particular,  as  an  — >  1,  these  have  the  limits 

Fo(an  — ►  1)  ->  0.707  +  0.680(1  -  an)  (3.10) 

Fi(an-^1)  ->  0.707  +  0.821(1  -  an) 

which  are  close  to  (3.6).  The  other  limit  for  an  produces 

A„(0;an  -»  0;  A)  -»  1.1215(3  -  A)  -  (8.446  -  5.806A)an.  (3.11) 

For  the  symmetric  double  crack  case,  Tada  et  al. ,  1985,  page  19.1  presents  Fn( 0;  an;  A) 
as  before  but  with 

F0(an)  =  0.5(3  -an)(l  +  1.243(1  -  an)3)  (3.12) 

Fi(an)  =  1  +  0.5(1  -an)  +  0.743(1  -  an)3. 

In  this  case,  as  an  — >  1,  both  F0  and  Fi  approach  the  same  limit  as  (3.8).  The 
an  — *  0  limit  is 

An(0;  an  ->  0;  A)  -*  1.1215(3  -  A)  -  (6.715  -  3.986A)cn.  (3.13) 

The  Fa  and  F\  functions  in  both  cases  are  accurate  to  within  1  %  according  to 
Tada  et  al. ,  1985. 

Tada  et  al.,  1985,  page  8.1  also  presents  the  edge  crack,  essentially  the  limit  of  the 
above  two  cases  as  an  — ■>  0,  without  the  stress  concentration  factor  of  (3  —  A).  For 
the  edge  crack,  A  has  no  effect,  and 

Ffse{ka  =  0)  =  1.1215.  (3.14) 

Later  arguments  for  small  ka  indicate  that,  for  all  geometries,  Fn  may  be  expanded 
in  the  series 

Fn(ka  — >  0;  an]  A)  — >  aF(an;  A)  -  bf(an ;  A )ka  +  cF(an ;  A )(ka)2,  (3.15) 

where  aF  gives  7^(0;  an;  A)  as  reproduced  from  Tada  et  al.,  1985  above.  A  value  may 
be  calculated  for  bF  in  the  edge  and  centre  crack  cases  following  the  self-consistent 
perturbation  approach  of  the  next  section.  For  the  hole  cases  it  is  calculated  numer¬ 
ically  by  fitting  (3.15)  to  Fn  calculated  for  small  ka.  Rose,  1987  gives  the  numerical 
value  cF  =  2.110  for  the  centre  crack  case. 

3.5  Self-Consistent  Perturbation  Theory  for  Small  ka 


There  are  three  ideas  underpinning  this  approach,  the  first  being  that  A'tip  may 
be  calculated  as  a  superposition  of  the  A't;p  values  produced  by  the  contributing 
stresses:  those  due  to  the  applied  load  and  the  springs.  The  second  is  that  those 
stresses  acting  closest  to  the  crack  tip  contribute  most  to  Kfip,  hence  the  crack 
opening  profile  needs  to  be  known  most  accurately  there  to  produce  an  accurate 
value  for  A\sf ,  the  contribution  of  the  springs.  Thirdly,  we  take  the  crack  profile 
to  be  the  same  as  the  unsprung  crack  with  the  same  loading  geometry,  but  scaled 
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(self-consistently)  to  give  the  correct  Kiiv  according  to  equation  (2.49).  These  ideas 
will  be  illustrated  for  the  centre  crack  case  which  can  be  treated  analytically. 

If  g{x\  a)  is  Ktip  for  a  crack  of  length  specified  by  a  (length  is  2 a  for  the  centre  crack) 
for  a  unit  applied  point  force  at  x  on  the  crack  face,  then  the  total  stress  intensity 
factor  is 

Ktip=  [  a(x)g(x;a)dx.  (3.16) 

J crack 

Here,  <r(x)  is  the  net  stress  on  the  crack  face  due  to  the  applied  load  and  the  springs. 
Using  the  first  idea  above,  it  may  be  split  as 

<r(x)  =  <r“(x)  -  i E'k6(x)  (3.17) 

where  equations  (2.3)  and  (2.27)  have  been  used.  This  leads  to 

Kt ip  =  Ko-^rE'k  f  S(x)g(x\a)dx  (3.18) 

2  J crack 

Ko  =  [  <7™y{x)g(x;a)dx. 

Jc  rack 

For  the  centre  crack,  cr™(x)  is  a  constant,  cr^. 


If  ka  is  small,  then  to  lowest  order  S(x)  will  be  the  same  as  if  there  were  no  springs. 
This  would  be  sufficient  for  it  enters  into  the  Ktip  correction  term  only.  Ideas  two 
and  three  above  suggest  that  we  can  do  better  than  this,  for  this  S(x)  is  known  to 
produce  K0  by  equation  (2.49).  For  the  centre  crack  and  uniaxial  remote  tension, 
Tada  et  al,  1985,  page  5.1a  gives 

4_oo  _  4/t°° 

6o(x)  =  — fr\ll  ~  (x/°)2  =  -~^r8\x)  (3-19) 

K0  =  CTyyV^a- 

We  replace  K0  in  the  second  equation  by  Ktip,  and  substitute  for  cr~  in  the  first, 
thereby  accomplishing  the  third  idea  above.  The  result  is 

Ktip  =  Ko - ~=-  /  6'(x)g(x]a)d(x/a).  (3.20) 

\/7ra  J  crack 


Substituting  in  S'(x)  from  above,  and 


(  x  2  1 

g(x\a)  =  —=—j^=^= 

1  -  ( x/a)2 


(3.21) 


from  Tada  et  al. ,  1985,  pages  5.11  and  5.11a  for  point  forces  applied  at  x  and  —x  in 
the  centre  crack,  the  equation  for  Ktip  reduces  to 


4  hn  /l 

K0  —  ——Ktip  I  d(x/a) 

7T  Jx/a=0 

K0(l  -  (4/tt )ka  +  (I6/7 r2)(ka)2  H - ). 


This  compares  to  the  Rose,  1987  form 


(3.22) 


Ktip/ Kq  =  1  -~ka  +  2.m(ka)2  + 

7 r 


(3.23) 


The  first  order  term  would  have  been  the  same  even  if  the  self-consistent  replacement 
of  Ko  by  Ktip  had  not  been  done,  but  that  would  not  have  given  a  second  order  term. 
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Figure  8:  Hole  co-ordinates  related  to  those  for  the  edge,  corresponding  to  an  — >  0. 
The  crack,  extending  from  R  =  b\  to  R  +  a  —  br  is  treated  like  an  edge  crack  from 
X  =  0  to  1.  Note  that  the  normalization  is  different  to  equation  (2.45)  but  that  it 
is  also  applied  to  the  double  crack  case. 


The  edge  crack  may  be  treated  similarly  as  the  8(x)  and  g(x ;  a)  functions  are  given 
(approximately)  by  Tada  et  al. ,  1985,  pages  8.1a,  8.3.  Just  presenting  the  equa¬ 
tions, 

8'(x)  =  a%J\  —  ( x/a)2D(x/a ) 

D(  X)  =  1.454  —  0.727A'  +  0.618A'2  —  0.224.Y3  (3.24) 

g(x;a)  =  "  , - 1==G(x/a ) 

V™/l  -  ( x/a )2 

G(  X)  =  1.3-0.3Ar5/4 

Aka  r1 

A'tip  =  K°-Yj2iteK'*J0  D{X)G(X)AX. 

The  integral  equals  1.4551  making  the  final  result 

Kt ip  =  cr^V7ra(l.l215  -  1.8516 ka  -f  3.0587(A:a)2  +  •••).  (3.25) 


The  other  cases  examined  here,  single  crack  from  a  hole  and  symmetric  double 
crack  from  a  hole,  can  not  be  treated  so  simply  because  the  functions  8(x)  and 
g(x;  a;  R)  are  not  known.  The  best  that  can  be  done  is  to  fit  appropriate  functions  to 
numerically  calculated  profiles  for  <f(:r),  and  to  seek  an  approximation  for  g(x:  a ;  R). 
Both  cracked  hole  cases  will  be  related  to  the  edge  crack  treatment  as  illustrated  in 
fig.  8.  This  is  the  an  — *  0  correspondence,  to  be  compared  to  fig.  7  which  was  for 
►  1-  the  case  of  uniaxial  tension,  A  =  0,  the  correspondence  may  be  summed 
up  as 


5(*)  = 


4<vg. 


1  -  A "D{X 


a  —  br  —  bi  —  (/?-)-  a)  —  R  (i 

X  =  (x  -  bi)/(br  -  bi)  =  (x  -  R)/a. 

Noting  the  original  normalizations  (2.45),  the  function  D( X)  is  obtained  using 

x  =  =  (^Tr) (A'  -  ^ where  Bi  = 1  <3 


(3.26) 


(X  —  1),  where  Bi  —  1 


(3.27) 
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Table  3:  The  cubic  function  D(X)  fitted  using  least  squares  to  the  numerical  crack 
profiles  for  various  an.  The  normalizations  for  X  and  D(X )  are  given  by  equations 
(3.27). 


an  D(X)  for  an  edge  crack 
edge  1.454  -  0.727Z  +  0.618X2  -  0.224V3 
3  x  edge  4.362  -  2T8lX  -f  1.854V2  -  0.672V3 
an  D(X )  for  a  single  crack  from  a  hole 
0.1  3.605  -  1.787V  +  1.399V2  -  0.499V3 
0.5  1.720  -  0.871V  +  0.701V2  -  0.245V3 
0.9  0.523  +  0.493V  -  0.313V2  +  0.083V3 
an  D(X)  for  a  double  crack  from  a  hole 
0.1  3.626  -  1.798V  +  1.408V2  -  0.502V3 
0.5  1.983  -  1.017V  +  0.767V2  -  0.262V3 
0.9  1.139  —  0.262V  +  0.330V2  —  0.153V3 
0.9*  1.137  -  0.257V  +  0.322V2  -  0.149V3~ 

*  Values  calculated  using  n  =  100  points  indicate  degree  of  convergence,  and  accu¬ 
racy  obtained  by  n  =  200. 


D(X) 


E'  R8n(X ) 
4  Ky  h*  ~  b‘ 


tn(X) 


Cubic  polynomials  were  then  fitted  to  D(X)  using  least  squares  as  listed  in  table  3. 
These  functions,  for  both  hole  cases,  are  plotted  in  fig.  9.  The  edge,  3  times  edge, 
and  centre  crack  functions  are  plotted  for  comparison.  The  functions  are  particularly 
smooth,  and  the  interpolations  should  be  comparable  in  accuracy  to  the  edge  case. 
Both  hole  cases  approach  3  times  the  edge  case  as  an  0,  but  only  the  double 
crack  case  approaches  the  centre  crack  at  the  other  limit.  This  is  because  the  single 
crack  case  approaches  a  centre  crack  of  full  length  a  rather  than  2 a  in  the  an  — »  1 
limit.  X  =  0  is  a  closed  end  of  the  crack  in  this  case  whereas  it  is  the  middle  of  the 
crack  in  the  double  crack  case.  For  the  single  crack  and  an  — *  0,  D(X  =  0)  — >  0 
and  D(X  =  1)  ->  l/y/2. 

The  crack  profile  may  also  be  expressed  using  the  correspondences  of  fig.  7  (c)  and 
(d).  For  the  single  crack,  this  gives  a  function  Z^.^)  which  is  almost  constant  over 
most  of  the  X2  range,  whilst  correspondence  7  (a)  diverged  for  X2  — >  0.  Still  on 
correspondence  7  (c),  equations  3.27  are  replaced  by 

—  2 x  -  a  _  2(1  —  an)  (  ^  2an  —  . 


d2{x2) 


=  — — — — —  (V  +  1)  -  1  =  -^-(X  -  1)  +  1 
2 R  +  a  2  —  an  v  ;  2-anv  ! 

E'  2 R6n(x)  1 

4<t-  2 R  +  a  s/l_x^ 


(3.28) 


These  functions  are  shown  in  fig.  10.  They  are  more  nearly  constant  than  D(X) 
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Figure  9:  The  function  D( X)  from  equation  (3.27)  for  an  =  0.1, 0.5,  and  0.9  com¬ 
pared  with  the  centre  and  edge  crack  functions.  Here  A  =  0  for  both  hole  cases. 


Figure  10:  The  functions  /^(A^),  equations  (3.28),  for  an  =  0.1,  0.5,  and  0.9  for  a 
single  crack  from  a  hole.  Note  the  expanded  D2  (vertical)  axis  in  each  case.  The 
arrows  indicate  X2  values  corresponding  to  X  =  0. 


(note  the  expanded  scale  on  fig.  10),  but  the  rapid  change  as  X2  — >  0  for  an  =  0.9 
may  cause  fitting  problems.  As  mentioned  earlier  though,  the  profile  has  to  be  most 
accurate  near  the  crack  tip.  In  this  case,  and  checks  would  have  to  be  made  to 
verify  it,  the  benefit  of  near-constancy  over  most  of  the  X2  range  may  override  the 
problems  at  X  —  0. 


Turning  to  the  double  crack  case,  the  preferred  correspondence  is  given  in  fig.  7  (d). 
In  this  case,  the  normalizations  are 


D3{  X3) 


—  =  (l-«n)*  =  flB(X-l)  +  l 
E'  R6n(x)  1 

4(Tw  r  +  a  yj\  -X3'2 

*  /1  _  X  sn(X3)  _  [Z  ,  //- 


(3.29) 


-«n) 


1 1  -  X* 


=  J(l+X)/(X-l  +  2/an)D(X). 


These  £>3(^3)  functions  are  plotted  in  fig.  11  for  the  same  an  as  before.  They  have 
no  advantages  over  the  functions  D(X)  which  are  thus  preferred  because  they  are 
well  behaved  as  an  — >  0  and  1. 


We  thus  have  D( X)  functions  for  both  hole  cases  that  are  as  reliable  as  D( X)  for  the 
edge  case,  but  we  still  need  functions  analogous  to  G(X )  for  the  edge  crack.  These 
are  not  so  easily  obtained,  as  Ktip  for  a  unit  point  force  applied  to  the  crack  faces  is 
required  here,  whereas  for  D(X),  the  information  was  already  available  through  the 
crack  profile.  A  separate  calculation  would  be  required,  beyond  the  scope  of  this 
piece  of  work.  We  may  expect  that  g(x]a;R)  for  the  single  crack  from  a  hole,  lies 
between  the  limiting  edge  and  centre  crack  cases  of  fig.  12.  These  correspondences 
have  been  presented  in  equations  (3.27)  and  (3.28),  so  only  the  normalizations  for 
the  g(x;a;R)  functions  need  be  given.  The  edge  case  was  presented  in  equations 
(3.24).  The  function  G(X)  given  there  should  be  an  upper  limit  for  the  hole  cases 
because  it  gives  Ktip  assuming  the  crack  is  from  an  edge.  With  either  hole  case,  the 
material  on  the  opposite  side  of  the  hole  to  Ktip  will  act  to  reduce  Artip  to  less  than 
this  estimate. 


For  both  hole  cases,  a  lower  bound  for  g(x;a;R )  is  found  by  assuming  the  hole  is 
part  of  a  centre  crack  as  in  fig.  7  (c)  and  (d).  This  is  because  removal  of  material  to 
form  the  hole  would  allow  some  relaxation  at  the  point  of  application  of  the  point 
force,  and  increase  I<tiP. 

Looking  at  the  single  crack  case  first,  the  approximation  requires  g(x]a]R )  for  a 
point  force  acting  on  a  centre  crack,  and  not  the  symmetric  pair  that  were  used  for 
the  centre  crack  case  (equation  3.21).  Tada  et  al,  1985,  page  5.10  gives,  for  a  centre 
crack  of  length  2 a  and  point  force  at  x, 

2  11 

gc{x-,a)  =  - -j=  =t(i  +  x/a).  (3.30) 

V™yf 1  -  {x/a)2  2V  > 

Here,  the  first  two  factors  are  the  function  from  (3.21),  separated  out  to  facilitate 
later  work.  We  use  the  same  substitutions  as  (3.28),  where  x/a  above  becomes  X2 
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Figure  11:  The  functions  D3(XS)  for  the  symmetric  double  crack  from  a  hole  case. 
The  functions  are  shown  for  an  =  0.1,  0.5  and  0.9.  These  are  not  really  any  more 
useful  than  the  D(X)  functions  of  fig.  9. 


ge(xe'>ae) 

a  e=a 
xe=x-R 


8c(Xc>‘ac) 

2R+a 

a°-  2 

xc=x-a/2 


Figure  12:  The  function  g(x\  a;  R)  for  the  single  crack  hole  case  would  be  expected 
to  lie  between  the  edge  and  centre  crack  limits.  ge(xe]  ae)  and  gc(xc ;  ac)  respectively, 
as  shown. 


and  a  is  replaced  by  (2 R  +  a)/ 2.  The  result  is 


g{x2)  = 


rx^{1+*2)- 


(3.31) 


To  enable  comparison  with  the  edge  crack  (maximum  bound),  the  edge-crack  nor¬ 
malization  is  introduced.  From  equation  3.24, 


g(x2)  = 


G2(X2 


:(l-X2) 


G2(X2) 


(3.32) 


I  a  1  -  X 

2(2/?  +  ci)  \  1  —  X2 


'(1+X)(1  +  A'2). 


The  relationship  between  X  and  X2  was  given  in  equations  (3.28).  Using  this, 
G2(X2)  can  be  expressed  as  a  function  of  X , 


G-2(X)  =  J~~^X  +  1)(X  -  2 +  2/an), 


(3.33) 


which  is  a  hyperbola.  It  is  plotted  for  several  an  along  with  G( X)  in  fig.  13. 
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X 


Figure  13:  The  functions  G(X)  (edge  approximation)  and  G2{X)  (approximating 
to  a  centre  crack)  for  the  single  crack  from  a  hole.  The  latter  function  is  shown 
for  several  an ,  and  is  part  of  a  hyperbola  in  each  case.  The  true  function  would  lie 
between  G(X)  and  G2(X)  for  each  value  of  an. 


It  is  clear  from  that  figure  that  these  bounds  are  too  wide  to  be  useful:  the  range  of 
possible  G  functions  is  so  great  that  first  order  coefficients  in  the  small  ka  expansion 
would  be  poorly  specified  here.  We  may  as  well  take  G(X)  =  1  in  the  integral 
replacing  (3.24)  for  the  single  crack  from  a  hole,  for  it  is  certainly  within  the  bounds. 
Given  this  uncertainty,  this  part  of  the  investigation  was  not  carried  any  further,  and 
G3(X()  functions  (see  equation  3.29)  for  the  double  crack  case  were  not  examined. 
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4  NUMERICAL  RESULTS  AND  INTERPOLATION  FORMULAE 

FOR  Ktip 

In  this  section,  the  numerical  results  for  Fn(ka ;  a„;  A)  will  be  presented.  They  have 
been  calculated  for  an  in  steps  of  0.1,  for  A  =  0  and  1,  and  for  both  hole  cases, 
(fully  bridged)  centre  and  edge  cracks.  These  are  shown  in  fig  14,  and  tabulated  in 
appendix  A.l.  In  many  cases,  it  may  be  sufficient  to  interpolate  in  this  table. 

A  simpler  method  of  presenting  and  using  this  data  is  to  construct  interpolation 
formulae  against  ka ,  given  an  and  A,  and  to  tabulate  only  the  parameters  of  these 
formulae.  Before  proceeding  to  these,  it  is  useful  to  consider  a  “map”  of  the  avail¬ 
ability  of  asymptotic  limits  against  which  the  numerical  results  may  be  verified. 
These  are  presented  in  fig  15.  They  also  serve  to  indicate  suitable  forms  for  the 
interpolation  formulae,  which  have  been  fitted  following  the  method  of  Rose,  1987. 
The  parameters  are  presented  in  fig  16  and  tabulated  in  appendix  A. 2. 

We  note,  as  did  Rose,  1987,  that 

tp  ( u„.  „  .  \  \ /  A)  bp(an,  A )ka  T  A)(^g)  ,  ka.  — *  0  ,  . 

J  -  1  *•(<.»;  JOA/K,  ka^oo  ^ 

so  that  an  appropriate  interpolation  formula  is 


F^(ka-an]  A)  = 


'  sp + ppka 
1  4-  qpka  +  rp(ka)2 ' 


where  the  (an;  A)  dependences  of  the  parameters  have  been  omitted  for  brevity.  On 
comparing  the  small  and  large  ka  expansions  of  this  with  the  previous  equation,  the 
four  parameters  pp  to  Sp  can  be  found  in  terms  of  ap  to  dp.  These  are 


sp  =  aF  rp 


Pf  =  d2Frp  qF 


2  spcp  —  Saptip 
2bpd2F  —  apsF 
2  aFbF  +  pF 


From  equation  (3.2), 


dF  -  (l  +  -(1  +  A)(l  -  an)2  +  -(1  -  A)(l  -  an)4^  (4.4) 

for  both  hole  geometries.  The  pp  to  sp  parameters  are  presented  in  fig.  16  for  both 
hole  cases  and  uniaxial  as  well  as  biaxial  loading. 


5  ASYMPTOTIC  LIMITS  FOR  CRACK  MOUTH  OPENINGS 


A  number  of  asymptotic  limits  may  be  examined  for  the  crack  mouth  openings, 
the  normalized  form  of  which  were  given  in  equation  (2.69).  These  are  similar  to 
those  for  Ktip,  but  not  so  extensive.  They  often  require  more  extreme  values  of  the 
parameters  before  becoming  “good”  approximations.  Some  difficulties,  in  addition 
to  those  noted  by  Rose,  1987,  were  experienced  in  fitting  interpolation  formulae  to 
Vn(ka;  an;  A).  These  functions  must  therefore  be  used  with  some  care. 
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Single  crack  from  a  hole  Double  crack  from  a  hole 


Figure  14:  Representative  plots  of  the  Fn{ka\an\  A)  functions  for  A  =  0  and  1,  for 
both  hole  cases  and  selected  values  of  an. 


_ ‘ _ 1(3.6) 

(3'13)  a=a/(a+R)  l3'8) 


Figure  15:  Diagram  showing  a  “map”  of  the  asymptotic  limits  to  Fn(ka.]  on;  A),  and 
which  equations  these  are. 


30 


0.0  0.2  0.4  0.6  0.8  1.0  0.0  0.2  0.4  0.6  0.8  1.0 

an=a/(a+R)  an=al(a+R) 


Figure  16:  Parameters  for  the  Fn  interpolation  formulae,  equations  (4.2). 


5.1  Stiff  Springs 


In  the  limit  of  stiff  springs,  away  from  the  crack  tip,  the  crack  opening  becomes  such 
that  “the  springs  carry  all  the  stress  which  would  have  existed  if  the  crack  were  not 
there.”  This  is  another  way  of  saying  that  Ktip  is  bounded  with  respect  to  increases 
in  a  by  the  reduction  effected  by  the  springs,  equation  (3.18). 

Assuming  large  kln  in  equation  (2.47),  the  second  term  dominates  over  the  first. 
Using  (2.1)  in  normalized  form,  we  obtain 

2t xkln8n{X)  =  -S%h(X).  (5.1) 

Vn(ka  — >  oo;an;  A)  is  found  for  X  —  Bi  using  (2.69),  the  second  equation  of  (2.47) 
and  (2.45),  for  both  hole  cases  as 

Vn(ka  — >•  oo;  an;  A)  -»  ^(3  -  A)-^.  (5.2) 

This  limit  requires  k  so  large  that  even  for  small  o,  ka  1.  If  a  is  too  small,  then 
the  criterion  stated  earlier,  “away  from  the  crack  tip”,  is  not  satisfied.  In  this  case 
the  material  just  beyond  the  crack  tip  will  be  carrying  some  of  the  stress  that  should 
be  carried  by  the  springs.  These  would  thus  not  be  extended  as  much  as  expected. 


5.2  Large  Hole  or  Short  Crack:  an  — >  0,  ka  =  0 


Both  hole  cases  tend,  in  the  an  — »  0  limit  with  no  springs,  to  an  edge  crack  with  an 
applied  stress  given  by  (3.3).  From  Tada  et  a/.,  1985,  page  8.1a, 

U„(0;  0;  A)  =  1.454(3  -  A),  ka  =  0.  (5.3) 


5.3  Small  R  or  Large  a:  an  — >  1  with  ka  =  0 


The  asymptote  in  the  limit  of  small  holes,  or  long  cracks,  is  not  easy  to  determine 
because  we  are  determining  behaviour  at  the  hole  boundary.  This  is  in  contrast 
to  the  Fn(ka;an  — »  1;  A)  limit  where  the  hole,  remote  from  the  crack  tip,  became 
increasingly  less  significant  as  an  — *  1. 

The  difficulty  is  best  exhibited  by  the  following  considerations.  Take  a  fixed  hole 
size  R,  and  a  series  of  increasingly  long  cracks  a.  For  stiff  springs  where  kR  >>  1 
and  ka  >>  1,  we  should  have  the  limit  (5.2).  In  the  case  of  no  springs,  as  even  weak 
springs  will  support  any  finite  stress  if  the  opening  is  wide  enough,  6(R)  would  be 
expected  to  increase  in  proportion  as  a  increases  for  the  double  crack  case.  For  the 
single  crack  case,  this  will  be  reduced  to  a  square  root  increase  by  the  material  on 
the  other  (uncracked)  side  of  the  hole.  An  approximation  for  both  hole  cases  is  to 
treat  the  hole  as  part  of  a  centre  crack  as  in  fig.  7  (c)  and  (d).  The  profile  of  an 
unsprung  centre  crack  was  given  in  equation  (3.19),  from  which  Vn(0;an  — >  1 ;  A)  is 
found  by  (2.69).  Care  is  required  with  the  change  of  variables  in  applying  (3.19)  to 
the  correspondences  of  fig.  7  (c)  and  (d). 
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Taking  the  single  crack  case  first,  a  becomes  (2 R  +  a)/ 2  and  x  i— >  x  —  a/2.  The 
opening  is  then  calculated  at  x  =  f?,  <5(1?)  =  4Sg?y  \j2aR.  This  leads,  as  an  — >  1,  to 

K(0;  an  ->  1;  A)  -  «  ^/2(1  -  a„)(i  +  (1  -  a„)/2  +  •■•).  (5.4) 

V 

For  the  double  crack,  we  use  a  i->  a  +  R  and  x  unchanged  but  set  equal  to  R  in 
(3.19),  whereupon 

14(0;  a„  -*  1;  A)  -»  ,/^S  a  1  +  (1  -  a„)  +  (1  -  a„)2/2  +  •  ■  ■ .  (5.5) 

V 


5.4  Large  a  with  Springs 


The  previous  section  examined  the  large  a  or  small  R  limit  in  the  absence  of  springs, 
indicating  different  behaviour  to  that  expected  when  finite  stiffness  springs  are 
present.  In  this  section,  the  integral  equation  (2.47)  is  re-examined  under  the  trans¬ 
formations  X  =  l/W  and  B  =  1/(7  as  Br  — >  oo.  In  this  case,  the  crack  opening  is 
no  longer  given  by  equation  (2.49),  but  needs  the  additional  term  M°°)>  the  crack 
opening  at  infinity: 

tw  An  S 

6n(W)  =  jCr  D(C )—  +  Moo),  Moo)  =  (5.6) 

In  this  equation,  CT  =  l/Br.  The  integral  equation  to  be  solved  then  becomes,  in 
the  infinite  crack  limit  (Cr  — »  0), 

/'  F(l/W,  C)^-iC  -  2 KkU  [  ^  =  -S„‘“(l/W),  (5.7) 

snuw)  =  I  Hw2 + ^4)  ■  A  = 0 

»  y  11  \  -SyW\  A  =  1 

It  is  not  necessary  to  invoke  the  second  change  of  variable  W  —  sin(t)  in  this  case, 
but  we  set  Q(C)  =  D(C)/C2,  and  the  crack  mouth  opening  becomes 

s^=^k+l'mia  (5'8> 

It  is  no  longer  possible  to  calculate  a  normalized  crack  mouth  opening  function  Vn 
so  comparisons  with  finite  cracks  will  be  made  using  Ml).  Similarly,  the  spring- 
stiffness  will  be  expressed  as  kR  rather  than  ka.  Apart  from  these  differences, 
numerical  solution  was  similar  to  that  described  previously.  The  results  are  pre¬ 
sented  in  appendix  A. 5  as  tables  of  Ml)  against  kR ,  calculated  from  above  for  the 
infinite  crack  cases,  and  from  (2.47)  for  finite  cracks.  An  interesting  feature  of  these 
results  is  that,  although  the  crack  mouth  opening  becomes  infinite  for  the  infinite 
crack  cases  as  kR  — >  0,  the  difference  Ml)  ~  M°°)  remains  finite. 

Figure  17  shows  how  the  crack  profile  near  the  crack  mouth  varies  as  the  crack  grows 
and  becomes  infinite  for  kR  =  1  in  the  single  crack,  A  =  0  case. 
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Figure  17:  Near-mouth  crack  profiles  for  a  single  crack  from  a  hole  in  uniaxial  (A  =  0) 
tension  for  varying  crack  lengths  and  kR  =  1.  (a)  to  (c)  show  the  dislocation  density 
and  (d)  to  (f)  the  near-mouth  crack  profile.  The  crack  tips  are  at  X  =  5,20,  and 
oo  for  (d),  (e)  and  (f)  respectively. 


5.5  Weak  Springs,  ka  — ►  0 


In  this  limit,  previous  arguments  for  Ktip  depended  on  the  assumption  that  the  crack 
profile  for  weak  springs  would  be  essentially  the  same  as  for  no  springs,  but  scaled 
down  to  produce  Ktw  instead  of  Kq,  the  stress  intensity  factor  for  no  springs.  In 
this  approximation,  the  crack  mouth  opening  8(R)  would  become  8(R)  =  60(R)-} :-R-. 
This  last  factor  was  expanded  as  a  series  in  ka,  evaluated  for  the  centre  and  edge 
cracks,  with  a  similar  expression  expected  for  the  hole  cases.  We  assume  therefore, 
that 

Vn (ka  — >  0;  an;  A)  — *  av(an:  A)  -  bv{an ;  A )ka  +  cv(an\  A )(ka)2.  (5.9) 

This  is  supported  by  noting  that  D(B)  in  equation  (2.47)  may  be  expanded  in 
powers  of  kln  for  ka  so  small  that  kR  is  too  (ln  =  R  for  holes,  a  for  edge  and  centre 
cracks) * 

D(B)  =  D(0\B)  -  klnD(1)(B)  +  (kln)2D{2)(B) - .  (5.10) 

Insertion  into  (2.47)  and  examination  of  the  powers  of  kln  produces 

(kin?  ■  fB' F(X,B)D<°\B)dB=-S~HX)  (5.11) 

J  Bi 

(kQ1  ■  [Br  F(X,  B)D™(B)dB  =  —2n  f*  D^(B)dB  =  -2t 

J  Bi  v  X 

The  first  of  these  is  just  the  unsprung  equation,  whilst  the  second  is  the  same 
equation  but  with  a  different  “applied  stress” .  It  would  thus  have  a  similarly  well- 
behaved  solution,  the  first  order  dislocation  density  D^(B).  Substituting  (5.10)  into 
the  second  equation  of  (2.49)  will  then  produce  a  corrected  crack  mouth  opening  of 
the  form  of  (5.9). 

6  NUMERICAL  RESULTS  AND  INTERPOLATION  FORMULAE 
FOR  Vn 


Numerical  values  for  Vn(ka\an\  A)  were  calculated  in  a  similar  manner  to  those  for 
Fn(ka\  an;  A).  Asymptote  (5.5)  was  excellent  for  A  =  1,  but  poor  for  A  =  0  as 

shown  in  table  4.  Equation  (5.4)  required  the  second  term  to  produce  a  satisfactory 

agreement  for  A  =  1,  but  was  also  poor  for  A  =  0. 

For  an  —>  0,  the  predicted  limit  by  (5.3)  is  Vn  =  2.908  for  A  =  1,  and  4.362  for  A  =  0. 

These  are  compared  with  calculated  values  in  table  5. 


Based  on  (5.2)  and  (5.9),  a  suitable  interpolation  function  for  Vn(ka]  an]  X)  was 
chosen  as 


K“(fca;a„;A)  = 

<Sv  =  ay 
pv  =  dyry 


sy  +  pyka  py(ka  —  Zq) 

1  +  qyka  +  ry(ka)2  ry(ka  —  z{)(ka  —  z2) 

rv  =  (by  —  SyCy)/(sy  —  bydy) 

qv  =  (pv  +  by)  /  Sy 
dy  =  (3  —  A)/2. 


(6.1) 


This  produced  good  interpolations  for  most  values  of  an  and  A  but,  unlike  equation 
(4.2)  for  Fn(ka]an;  A),  the  denominator  here  had  a  root  in  the  ka  >  0  range  for 
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Table  4:  Testing  the  Vn(0;an  — >  1;  A)  asymptotes,  equations  (5.4)  and  (5.5). 


single 

crack  from 

a  hole 

&n 

14(0;  an;0) 

K(0  ;an;  1) 

1  term 

2  term 

1.1045 

0.9502 

■SBa 

0.891 

BIS 

0.8244 

0.7354 

0.696 

0.5333 

0.4954 

0.470 

0.3605 

0.3434 

0.324 

Bill 

0.2220 

0.2164 

0.202 

double  crack  from 

a  hole 

an 

K(0;  an;  0) 

14(0  ;a„;  1) 

linear 

quadratic 

EH 

1.5138 

1.3286 

i 

1.3222 

1.2108 

BB 

1.1516 

1.1021 

1.100 

1.105 

0.95 

1.0728 

1.0504 

1.050 

1.051 

0.98 

1.0267 

1.0200 

1.020 

1.020 

Table  5:  Testing  the  14(0;  an  —*  0;  A)  asymptote,  equation  (5.3). 


single  crack  from  a  hole 

double  crack  from  a  hole 

an 

K(0;a„;0) 

14(0;  an;  1) 

14(0;  an;0) 

K(0;a„;  1) 

0.1 

3.6077 

2.5029 

3.6292 

2.5180 

0.05 

3.9559 

2.6909 

3.9619 

2.6951 

0.02 

4.1834 

2.8117 

4.1845 

2.8124 

asymptote 

4.362 

4.362 

Table  6:  Testing  the  Vn(ka  — >  oo;an;  A)  asymptote,  equation  (5.2).  A  few  values  of 
an  have  been  included  to  show  that  the  limit  requires  higher  ka  as  an  increases. 


double  crack  from  a  hole 

ka 

Vn(ka ;  0.1;  0) 

Vn(k :a;  0.5;  0) 

Vn(ka- 0.9;  0) 

Vn (ka:  0.1;  1) 

■ 

1.1196 

0.7826 

0.4681 

0.7673 

0.4560 

0.3652 

0.2318 

0.3093 

0.15 

0.10 

0.1463 

0.1316 

0.0944 

0.0982 

0.05 

0.033 

0.0496 

0.0474 

0.0386 

0.0332 
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some  an.  This  is  shown  in  fig.  18(b)  for  the  single  crack  from  a  hole  case.  In 
addition  to  this,  the  pv,qv  and  ry  parameters  diverged  as  shown  in  fig.  18(a)  at  an 
an  dependent  on  the  hole  case  and  A. 

In  order  to  overcome  these  difficulties,  two  other  interpolation  formulae  were  exam¬ 
ined.  The  first,  more  approximate  one,  was  suggested  by  the  near-coincidence  of  zq 
and  z2  (equation  (6.1))  in  the  a„  range  causing  difficulty.  This  function  contains 
only  two  parameters, 


V:™(ka]an]  A)  = 


ka  T  (dy /ay) 


ay ,  ka  — >  0 
dy/ka,  ka  — >  oo 


which  correctly  matches  the  ka  =  0  and  ka  — >  oo  limits. 


An  improvement  is  the  three  parameter  “alternativee”  function 


V:x\ka-an-  A)  =  J- 


•Sa  _  -Sa/^a 

+  qaka  +  ra(ka )2  \  (ka  -  r-y)(ka  -  r2 ) 


sa  =  ay  ra  =  ( ay/dy )2  (6.3) 

qa  =  2 by /ay  ca  =  av(3q2a/4  -  ra)/2 

which  also  matches  the  gradient  as  ka  — >  0.  It  leads  to  the  coefficient 

ca  for  the  (ka)2  term,  which  may  be  compared  to  the  true  (numerically  obtained) 
value  cy. 


Given  that  by  is  obtained  numerically,  a  preferred  method  of  obtaining  the  three 
parameters  was  to  fit  to  ay,  dy  and  the  value  of  Vn  at  an  intermediate  value  of 
ka  =  ka0.  The  parameters  for  this  “new”  function  were: 

sn-a\  rn-  (ay/dy)2  (6.4) 

1  «sn 

Vn  —  Oi  ^ri)  —  T  T  (fcflo)  1 

kcio  \_vn 

The  parameters  sn,rn,  and  qn  for  A  =  0  and  1  for  both  hole  cases  are  shown  in 
fig.  19.  These  parameters  vary  smoothly  with  an,  in  contrast  to  those  for  V(Jnt.  The 
(denominator)  roots  rx  and  r2  also  cause  no  problem  because  they  occur  (when  real) 
at  negative  ka  only. 


Comparisons  between  the  various  interpolation  functions  and  the  calculated  values 
of  Vn(ka;  an\  A)  are  presented  in  fig.  20  for  the  single  crack  case  with  an  =  0.9,  A  =  0. 

is  excellent  for  ka  <  2  whilst  Vnapp  is  in  error  by  about  25%;  too  much  to  be 
useful.  The  situation  is  similar  for  an  =  0.5,  shown  in  fig.  21.  The  errors  for  all 
functions  are  reduced  by  a  factor  of  about  5,  and  the  alternatives  to  Fnmt  are  close 
enough  to  be  considered  useful. 


When  an  —  0.3,  the  Vnmt  function  has  a  singularity  in  it,  indicated  in  fig.  22  (a).  It 
causes  a  divergence  of  the  error  for  this  function,  shown  in  fig.  22  (b).  In  this  case, 
the  Vnalt  and  I4new  functions  are  actually  more  accurate  than  V),mt.  V(Jnt  possesses 
this  singularity  for  an  in  the  approximate  range  0.26  <  an  <  0.37. 
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0.0  0.2  0.4  0.6  0.8  1.0  0.0  0.2  0.4  0.6  0.8 

an-a/(a+R)  an=a/(a+R) 

Figure  19:  Plot  of  the  parameters  for  V^ew,  equation  (6.4),  against  an  with  A 
for  both  hole  cases. 


V  for  single  crack,  an-0.5,  h-0 


AJ(ka) 


(a) 


parabola\ 


n  n 


Percent  error  for  interpolations 


parabola/ 


d/(ka)\  (b) 


Figure  21:  Comparison  of  the  accuracy  of  the  various  Vn  interpolating  functions  for 
an  =  0.5.  Note  the  overall  reduction  in  error  compared  to  fig.  20. 


Vn  for  single  crack,  an-0.3,  X-0 


Figure  23:  The  symmetric  double  crack  case  (a)  compared  to  the  partially  bridged 
centre  crack  (b)  of  Rose,  1987.  The  variables  in  (b)  are  those  used  by  Rose,  1987. 

7  COMPARISON  OF  Ktip  FOR  THE  SYMMETRICALLY  CRACKED 
HOLE  AND  THE  PARTIALLY-BRIDGED  CENTRE  CRACK 


An  interesting  comparison  exists  between  Kt ;p  for  the  symmetrically  (double)  cracked 
hole  and  partially  bridged  centre  crack  cases.  In  particular,  the  question  of  which 
has  a  lower  I\t jp  arises,  and  thus  whether  drilling  out  the  unbridged  portion  of  the 
centre  crack  may  be  advantageous. 

Figure  23  shows  the  variables  used  here  (equation  2.68),  and  those  used  by  Rose,  1987 
for  the  partially  bridged  crack,  his  equations  15  (a)  and  (b).  Care  is  required  in  the 
comparison,  for  Rose,  1987  normalizes  Ktip  based  on  a  in  fig.  24  (b),  equivalent  to 
R  +  a  in  (24)  (a).  This  means  that  the  Fn(ka\  an\  0)  function  of  equation  (2.68) 

must  be  multiplied  by  a/(R  +  a )  =  ^/a. ~  before  comparison.  In  addition,  c/a  is 

equivalent  to  1  —  an. 

For  no  springs,  Rose,  1987  has  the  function  F(0,c/a)  =  1,  whilst  the  function  here 
will  be  v/a)fF'n(0;  an;  0)  =  y/dFaap{an\  0)  from  equation  4.1.  These  are  compared  in 
fig.  24  (a).  The  hole  case  has  a  lower  Ktip  for  c/a  >  0.82,  or  an  <  0.18. 

For  stiff  springs,  Rose,  1987,  equations  58  (b)  and  24  indicate  that  \[C  from 

F{kl-c/a)  =  y/C/Vkl,  y/C=  -  (7  1) 

+  c/a 

should  be  compared  with  ^/a^dp  from  (4.1).  This  is  done  in  fig.  24  (b)  where,  again 
as  c/a  — ►  1,  the  hole  case  becomes  more  favourable  with  a  lower  Ktip.  It  is  thus 
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better  (if  lowering  Kt\P  is  desirable)  to  drill  out  the  unbridged  portion  when  a  crack 
is  bridged  near  the  tips  only. 

In  fig.  25,  the  value  of  an  which  makes  Ktjp  equal  in  the  two  cases  is  plotted  as  a 
function  of  the  spring  stiffness  ka  ( kl ). 

8  CONCLUSIONS 


Values  of  normalized  crack  tip  stress  intensity  factor  from  equation  (2.68),  and 
normalized  crack  mouth  opening  (2.69),  along  the  lines  of  Rose,  1987  have  been 
presented.  These  are  for  the  cases  of  a  single  crack  or  symmetric  double  cracks 
emanating  from  a  hole  in  an  infinite  plate.  Both  uniaxial  and  biaxial  uniform  remote 
tensions  were  treated. 

Interpolation  functions  with  respect  to  spring  stiffness  have  been  provided  wherein 
the  function  parameters  depend  on  the  crack  length  relative  to  hole  radius,  an,  and 
the  biaxiality  A  of  the  remote  loading. 

Difficulties  were  experienced  with  the  interpolations  for  crack  mouth  opening  for 
some  values  of  an.  In  these  cases,  an  alternative  function  avoided  the  problem,  but 
this  function  was  less  accurate  away  from  these  an  values. 
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No  springs 


Rose(87) 

5xdifference 


a  -1-c/a 

n 


Figure  24:  Comparison  of  normalized  Ktjp  values  for  the  symmetric  double  crack 
hole  case,  and  the  partially  bridged  centre  crack,  (a)  is  for  no  springs  and  (b)  for 
stiff  springs. 


ka 


Figure  25:  Values  of  an  —  1  —  c/a  which  produce  the  same  Kttp  for  the  double  crack 
hole  and  partially  bridged  centre  crack  cases,  as  a  function  of  the  spring  stiffness. 
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A  TABLES  OF  VALUES 

Here  are  tabulated  the  various  numerical  results  which  are  provided  for  completeness 
and  reference  purposes.  The  tables  may  be  used  as  a  more  direct  but  less  accurate 
method  of  obtaining  Fn(^ka}  an)  A)  and  Vn[ka',  on;  A)  by  looking  up  the  values  for  an 
and  ka  closest  to  those  needed. 

A.l  Fn{ka\  an]  A) 

All  values  of  Fn(^ka ,  A)  presented  here  are  calculated  using  n  —  200  points  except 

for  the  few  values  for  case  (1,0)  indicated  by  #  n=100. 


#  fnkan  cases  here 

#  case=  edge ,  3*ed 

an\ka  .01  .1 


an\ka 

.05 

.10 

.20 

.40 

.60 

.80 

.90 

.95 

#  n=  100 
an\ka 
.10 
.40 
.60 
.90 


idge ,  3*edge, 
.01  .02 

centre 

.05 

.10 

.20 

.50 

.1029 

1 .0856 

1 . 0377 

.9690 

.8622 

.6705 

' .  3087 

3.2568 

3.1131 

2  .9070 

2.5866 

2.0115 

.9875 

.9754 

.9412 

.8907 

.8083 

.6483 

fd 

©  o 

.02 

.05 

.10 

.20 

.50 

.9732 

2.9283 

2.8033 

2.6231 

2.3399 

1.8241 

.6785 

2.6393 

2.6301 

2.3717 

2.1207 

1.6571 

.1901 

2 . 1600 

2.0757 

1 .9522 

1.7532 

1 . 3767 

.6177 

1 . 4993 

1.4472 

1.3696 

1.2413 

.9875 

.1150 

1.1034 

1.0702 

1.0199 

.9350 

.7598 

.8699 

.8624 

.8411 

.8084 

.7516 

.6292 

.7810 

.7751 

.7583 

.7322 

.6864 

.5847 

.7412 

.7361 

.7213 

.6983 

.6575 

.5654 

.01 

.02 

.05 

.10 

.20 

.50 

2.6769  2.6378  2.5289  2.3708  2.1201  1.6570 

1.5171  1.4987  1.4467  1.3692  1.2410  .9875 

1.1147  1.1031  1.0699  1.0197  .9348  .7598 

.7810  .7752  .7583  .7323  .6865  .5847 


i  case=  (1,1) 


4987 

1.4467  1.3692  1 

.2410 

.9875 

1031 

1.0699  1.0197 

.9348 

.7598 

7752 

.7583  .7323 

.6865 

.5847 

.02 

.05  .10 

.20 

.50 

2.0645  2.0337  1.9479  1 
1.9362  1.9085  1.8314  1. 
1.7110  1.6886  1.6255  1 


8243  1 .6298  1.2749 
7195  1.5419  1.2127 
5329  1.3834  1.0986 


.40 

1.3553  1.3399 

1.2962  1 

.2309 

1.1226 

.9064 

.60 

1.0863  1.0755 

1.0446 

.9979 

.9187 

.7542 

.80 

.8747  .8673 

.8462 

.8137 

.7574 

.6355 

.90 

.7846  .7788 

.7619 

.7357 

.6898 

.5877 

.95 

#  case= 

.7428  .7377 

(2,0) 

.7228 

.6998 

.6589 

.5666 

an\ka 

.01  .02 

.05 

.10 

.20 

.50 

2.9777  2.9326  2.8072  2.6264  2.3424  1.8255  1.4053 
2.6942  2.6545  2.5439  2.3836  2.1297  1.6621  1.2778 
2.2394  2.2078  2.1193  1.9899  1.7825  1.3930  1.0675 
1.6444  1.6223  1.5602  1.4684  1.3190  1.0317  .7866 
1.3054  1.2885  1.2408  1.1699  1.0537  .8276  .6325 
1.10801.09411.0549  .9967  .9011  .7143  .5519 
1.04111.0282  .9919  .9380  .8497  .6774  .5273 
1.01281.0004  .9652  .9132  .8279  .6620  .5173 


.60 

1.3054 

1.2885 

1.2408  1 

.1699  1 

.0537 

.8276 

.80 

1.1080 

1.0941 

1.0549 

.9967 

.9011 

.7143 

.90 

1.0411 

1.0282 

.9919 

.9380 

.8497 

.6774 

.95 

1.0128 

1.0004 

.9652 

.9132 

.8279 

.6620 

#  case= 
an\ka 

(2,1) 

.01 

.02 

.05 

.10 

.20 

.50 

2.0676  2.0366  1.9506  1 . 
1.9472  1.9192  1.8411  1 . 
1.7477  1.7240  1.6578  1 . 


1.6315  1.2758 
1.5482  1.2162 
1.4051  1.1108 


1.4590 

1.4405 

1.3886  1.3117 

1.1862 

.9425 

1.2561 

1.2405 

1.1967  1.1315 

1.0244 

.8143 

1.1028 

1.0892 

1.0508  .9939 

.9002 

.7166 

1.0408 

1.0280 

.9918  .9382 

.8504 

.6787 

1.0131 

1.0006 

.9655  .9135 

.8284 

.6625 

2.00 

5.00 

10.00 

20.00 

50.00 

.3821 

.2482 

.1771 

.1258 

.0799 

.1463 

.7446 

.5313 

.3774 

.2397 

.3801 

.2480 

.1770 

.1258 

.0799 

2.00 

5.00 

10.00 

20.00 

50.00 

.0356 

.6693 

.4760 

.3374 

.2138 

.9376 

.6030 

.4276 

.3025 

.1914 

.7751 

.4945 

.3490 

.2461 

.1553 

.5583 

.3542 

.2490 

.1752 

.1105 

.4427 

.2840 

.2008 

.1418 

.0897 

.3863 

.2538 

.1812 

.1287 

.0817 

.3706 

.2467 

.1772 

.1262 

.0802 

.3647 

.2445 

.1760 

.1256 

.0799 

2.00 

5.00 

10.00 

20.00 

50.00 

.9379 

.6035 

.4282 

.3032 

.1923 

.5585 

.3544 

.2494 

.1757 

.1110 

.4428 

.2843 

.2011 

.1422 

.0901 

.3708 

.2470 

.1775 

.1266 

.0807 

2.00 

5.00 

10.00 

20.00 

50.00 

.7294 

.4735 

.3375 

.2396 

.1520 

.6963 

.4516 

.3217 

.2282 

.1447 

.6348 

.4112 

.2924 

.2072 

.1313 

.5303 

.3430 

.2434 

.1722 

.1090 

.4495 

.2918 

.2073 

.1468 

.0929 

3921 

.2580 

.  1844 

.1310 

.0831 

3726 

.2480 

.1781 

.1268 

.0806 

3653 

.2449 

.1763 

.1257 

.0800 

2.00 

5.00 

10.00 

20.00 

50.00 

0360 

.6694 

.4760 

.3374 

.2138 

9388 

.6034 

.4277 

.3025 

.1914 

7789 

.4956 

.3494 

.2462 

.1554 

5684 

.3570 

.2501 

.1756 

.1106 

4577 

.2882 

.2023 

.1424 

.0898 

4045 

.2586 

.1829 

.1293 

.0819 

3899 

.2517 

.1789 

.1268 

.0804 

3844 

.2495 

.  1777 

.1261 

.0800 

.00 

5.00 

10.00 

20.00 

50.00 

7296 

.4735 

.3375 

.2396 

.1520 

6971 

.4519 

.3218 

.2283 

.1447 

6376 

.4120 

.2927 

.2073 

.1313 

5385 

.3453 

.2443 

.1725 

.1091 

4627 

.2954 

.2086 

.1472 

.0930 

4093 

.2626 

.1860 

.1315 

.0833 

3914 

.2529 

.1797 

.1274 

.0808 

3849 

.2498 

.1779 

.1263 

.0801 
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A. 2  Parameters  for  Fn  interpolations 


Below  are  tabulated  the  parameters  for  Fn(ka\  an\  A)  from  equation  4.2.  These  are  pp 
to  sp  together  with  the  roots  zi,  z2  of  the  denominator  and  zero  z0  of  the  numerator. 
Complex  roots  are  indicated  by  z\  —  z2  —  0. 


#  Fnint  cases  here 

#  case=  (1,0) 


# 

an 

s 

r 

P 

q 

zl 

z2 

zero 

.05 

9.120823 

5.048863 

11.482718 

4.423212 

.0000 

.0000 

-.7943 

.10 

7.393766 

4.228707 

7.683244 

4.090889 

.0000 

.0000 

-.9623 

.20 

4.933661 

3.320904 

3.955478 

3.654874 

-.5092 

-.5914 

-1.2473 

.21 

4.743758 

3.259147 

3.730515 

3.621078 

-.5133 

-.5977 

-1.2716 

to 

to 

4.562223 

3.200843 

3.522672 

3.588530 

-.5180 

-.6031 

-1.2951 

.23 

4.388668 

3.145690 

3.330386 

3.557139 

-.5232 

-.6076 

-1.3178 

.24 

4.222724 

3.093417 

3.152261 

3.526824 

-.5289 

-.6112 

-1.3396 

.25 

4.064039 

3.043781 

2.987054 

3.497509 

-.5352 

-.6139 

-1.3606 

.26 

3.912280 

2.996563 

2.833648 

3.469125 

-.5421 

-.6156 

-1.3807 

.27 

3.767128 

2.951563 

2.691041 

3.441608 

-.5498 

-.6162 

-1.3999 

.28 

3.628280 

2.908600 

2.558333 

3.414902 

-.5587 

-.6154 

-1.4182 

.29 

3.495447 

2.867509 

2.434709 

3.388953 

-.5696 

-.6123 

-1.4357 

.30 

3.368355 

2.828141 

2.319436 

3.363711 

-.5868 

-.6026 

-1.4522 

.31 

3 . 246740 

2.790357 

2.211847 

3.339131 

.0000 

.0000 

-1.4679 

.32 

3.130354 

2.754032 

2. 111340 

3.315170 

.0000 

.0000 

-1.4826 

.33 

3.018958 

2.719050 

2.017364 

3.291789 

.0000 

.0000 

-1.4965 

.34 

2.912326 

2.685305 

1.929421 

3.268952 

.0000 

.0000 

-1.5094 

.35 

2.810242 

2.652698 

1.847055 

3.246622 

.0000 

.0000 

-1.5215 

.36 

2.712499 

2.621139 

1.769849 

3.224769 

.0000 

.0000 

-1.5326 

.37 

2.618900 

2.590544 

1.697423 

3.203361 

.0000 

.0000 

-1.5429 

.38 

2.529259 

2.560834 

1.629428 

3.182371 

.0000 

.0000 

-1.5522 

.39 

2.443397 

2.531938 

1.565544 

3.161770 

.0000 

.0000 

-1.5607 

.40 

2.361143 

2.503788 

1.505477 

3.141533 

.0000 

.0000 

-1.5684 

.50 

1.704741 

2.251269 

1 . 064405 

2 . 954300 

.0000 

.0000 

-1.6016 

.60 

1.270007 

2.023035 

.805467 

2.781570 

.0000 

.0000 

-1.5767 

.70 

.975656 

1.789572 

.636609 

2.606221 

.0000 

.0000 

-1.5326 

CO 

o 

.769922 

1.527659 

.508298 

2.408181 

.0000 

.0000 

-1.5147 

.90 

.619217 

1.219432 

.392165 

2.159898 

.0000 

.0000 

-1.5790 

.95 

.557087 

1.046475 

.333943 

2.005735 

.0000 

.0000 

-1.6682 

# 

n=  100  below 

# 

an 

s 

r 

P 

q 

zl 

z2 

zero 

.10 

7.384899 

4.209123 

7.647661 

4.082275 

.0000 

.0000 

-.9656 

.50 

1.703599 

2.245012 

1.061446 

2.951274 

.0000 

.0000 

-1.6050 

.90 

.619314 

1.215541 

.390914 

2.157939 

.0000 

.0000 

-1.5843 
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# 

case= 

(1,1) 

# 

an 

s 

r 

P 

.05 

4.395849 

5.406093 

6.228493 

.10 

3.860911 

4.616570 

4.814229 

.20 

3.007906 

3.541231 

3.031741 

.30 

2.368473 

2.874534 

2.031375 

.40 

1.880227 

2.434853 

1.433510 

.50 

1.501610 

2.123572 

1.056178 

.60 

1.204138 

1.882849 

.806458 

.70 

.967810 

1.673547 

.632909 

.80 

.778273 

1.461311 

.503106 

.90 

.625040 

1.204982 

.391267 

.95 

.559528 

1.045490 

.334456 

# 

case= 

(2,0) 

# 

an 

s 

r 

P 

.05 

9.148752 

5.061794 

11.512127 

.10 

7.482468 

4.277533 

7.771956 

.20 

5.162956 

3.482781 

4.148287 

CO 

o 

3.711968 

3.135514 

2.571521 

.40 

2.779610 

2.980795 

1.792291 

.50 

2.165132 

2.921711 

1.381391 

.60 

1.750041 

2.915672 

1.160869 

.70 

1.462874 

2.946630 

1.048212 

.80 

1.259535 

3.019998 

1.004844 

.90 

1.111635 

3.173657 

1.020638 

.95 

1.052047 

3.309617 

1.056139 

# 

case= 

(2,1) 

# 

an 

s 

r 

P 

.05 

4.409128 

5.418157 

6.242393 

.10 

3.905996 

4.663816 

4.863497 

.20 

3.140531 

3.706863 

3.173543 

.30 

2.692679 

3.192142 

2.255822 

.40 

2.184530 

2.922182 

1.720423 

.50 

1.869219 

2.797765 

1.391494 

.60 

1.618146 

2.768740 

1.185901 

.70 

1.414148 

2.812244 

1.063545 

.80 

1.247149 

2.926437 

1.007525 

.90 

1.110828 

3.136583 

1.018473 

.95 

1.052559 

3.297806 

1.054980 

q 

zl 

z2 

zero 

4.542905 

.0000 

.0000 

-.7058 

4.227434 

.0000 

.0000 

-.8020 

3.740225 

.0000 

.0000 

-.9921 

3.384819 

.0000 

.0000 

-1.1659 

3.113904 

.0000 

.0000 

-1.3116 

2.897284 

.0000 

.0000 

-1.4217 

2.713650 

.0000 

.0000 

-1.4931 

2.544725 

.0000 

.0000 

-1.5291 

2.368870 

.0000 

.0000 

-1.5469 

2.149893 

.0000 

.0000 

-1.5975 

2.004897 

.0000 

.0000 

-1.6729 

q 

zl 

z2 

zero 

4.428970 

.0000 

.0000 

-.7947 

4.114280 

.0000 

.0000 

-.9628 

3.741105 

-.5006 

-.5736 

-1.2446 

3.539054 

.0000 

.0000 

-1.4435 

3.426222 

.0000 

.0000 

-1.5509 

3.367258 

.0000 

.0000 

-1.5674 

3.344799 

.0000 

.0000 

-1.5075 

3.351640 

.0000 

.0000 

-1.3956 

3.390039 

.0000 

.0000 

-1.2535 

3.477592 

.0000 

.0000 

-1.0892 

3.553635 

.0000 

.0000 

-.9961 

q 

zl 

z2 

zero 

4.548142 

.0000 

.0000 

-.7063 

4.249244 

.0000 

.0000 

-.8031 

3.825469 

.0000 

.0000 

-.9896 

3.563108 

.0000 

.0000 

-1.1493 

3.406087 

.0000 

.0000 

-1.2698 

3.321233 

.0000 

.0000 

-1.3433 

3.289304 

.0000 

.0000 

-1.3645 

3.300560 

.0000 

.0000 

-1.3297 

3.354364 

.0000 

.0000 

-1.2378 

3.463413 

.0000 

.0000 

-1.0907 

3.549074 

.0000 

.0000 

-.9977 
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A. 3 


Vn (&<z,  ctn,  A) 


Below  are  presented  the  numerically  obtained  values  of  Vn(ka\  an;  A). 


#  vnkan  cases  here 

*  case=  edge,  3»edge,  centre 


an\ka 

.01 

.02 

.05 

.10 

.20 

.50 

1.00 

2.00 

5.00 

10.00 

20.00 

50.00 

edge 

1.4153 

1.3830 

1 . 2941 

1.1680 

.9753 

.6462 

.4077 

.2306 

.0981 

.0497 

.0250 

.0100 

3* edge 

4,2459 

4.1490 

3 . 8823 

3 . 5040 

2.9259 

1.9386 

1.2231 

.6918 

.2943 

.1491 

.0750 

.0300 

centre 

.9823 

.9652 

.9171 

.8466 

.7331 

.5205 

.3485 

.2081 

.0932 

.0483 

.0246 

.0099 

#  case5 
an\ka 

(1,0) 

.01 

.02 

.05 

.10 

.20 

.50 

1.00 

2.00 

5.00 

10.00 

20.00 

50.00 

.05 

3.8688 

3.7853 

3.5543 

3.2243 

2.7147 

1 . 8273 

1.1692 

.6704 

.2894 

.1477 

.0745 

.0299 

.10 

3.5325 

3.4602 

3.2597 

2.9712 

2.5210 

1.7227 

1.1174 

.6494 

.2844 

.1462 

.0741 

.0298 

.20 

2.9519 

2.8976 

2.7459 

2.5249 

2.1733 

1.5284 

1.0182 

.6078 

.2740 

.1430 

.0731 

.0297 

.40 

2.0559 

2.0252 

1 . 9386 

1 . 8097 

1 . 5978 

1 . 1846 

.8306 

.5230 

.2507 

.1353 

.0708 

.0292 

.60 

1.3818 

1.3653 

1.3179 

1 . 2463 

1.1252 

.8757 

.6456 

.4297 

.2208 

.1243 

.0671 

.0285 

.80 

.8165 

.8089 

.7867 

.7527 

.6938 

.5660 

.4395 

.3113 

.1748 

.1049 

.0598 

.0267 

.90 

.5289 

.5246 

.5123 

.4931 

.4593 

.3843 

.3072 

.2260 

.1349 

.0853 

.0512 

.0243 

.95 

.3578 

.3551 

.3475 

.3356 

.3146 

.2671 

.2172 

.1635 

.1016 

.0669 

.0420 

.0212 

*  n=  100 
an\ka  .01 

.02 

.05 

.10 

.20 

.50 

1.00 

2.00 

5.00 

10.00 

20.00 

50.00 

.10 

3.5212 

3.4493 

3 . 2498 

2.9627 

2.5145 

1.7192 

1.1157 

.6487 

.2842 

.1462 

.0740 

.0298 

.40 

2.0504 

2.0199 

1 . 9335 

1.8051 

1.5939 

1.1818 

.8288 

.5220 

.2503 

.1351 

.0707 

.0292 

.60 

1.3789 

1.3623 

1.3151 

1.2435 

1.1226 

.8736 

.6439 

.4285 

.2202 

.1241 

.0670 

.0284 

.90 

.5293 

.5250 

.5125 

.4932 

.4593 

.3840 

.3066 

.2253 

.1343 

.0848 

.0509 

.0241 

#  case5 
an\ka 

(1,1) 

.01 

.02 

.05 

.10 

.20 

.50 

1.00 

2.00 

5.00 

10.00 

20.00 

50.00 

.05 

2.6314 

2.5744 

2.4167 

2.1914 

1.8435 

1 . 2384 

.7903 

.4517 

.1941 

.0988 

.0498 

.0200 

.10 

2.4503 

2.3997 

2.2595 

2.0578 

1.7432 

1.1864 

.7657 

.4422 

.1920 

.0982 

.0496 

.0199 

.20 

2.1268 

2.0870 

1.9758 

1.8139 

1.5567 

1.0864 

.7168 

.4226 

.1874 

.0969 

.0492 

.0199 

.40 

1.5921 

1 . 5674 

1 . 4977 

1.3942 

1.2245 

.8956 

.6174 

.3801 

.1766 

.0935 

.0482 

.0197 

.60 

1.1471 

1.1325 

1 . 0908 

1.0279 

.9219 

.7053 

.5087 

.3288 

.1617 

.0884 

.0467 

.0194 

.80 

.7279 

.7206 

.6994 

.6670 

.6109 

.4903 

.3724 

.2557 

.1363 

.0786 

.0432 

.0186 

.90 

.4911 

.4869 

.4747 

.4557 

.4226 

.3490 

.2740 

.1961 

.1114 

.0675 

.0388 

.0175 

.95 

.3408 

.3381 

.3305 

.3187 

.2977 

.2504 

.2010 

.1481 

.0884 

.0558 

.0336 

.0160 

*  case5 
an\ka 

(2,0) 

.01 

.02 

.05 

.10 

.20 

.50 

1.00 

2.00 

5.00 

10.00 

20.00 

50.00 

.05 

3 . 8745 

3.7908 

3.5592 

3.2284 

2.7176 

1.8286 

1 . 1697 

.6706 

.2894 

.1477 

.0745 

.0299 

.10 

3.5531 

3.4800 

3 . 2773 

2.9859 

2.5317 

1.7278 

1.1196 

.6502 

.2846 

.1463 

.0741 

.0299 

.20 

3.0188 

2.9621 

2 . 8038 

2.5740 

2 . 2099 

1.5466 

1.0264 

.6107 

.2747 

.1432 

.0732 

.0297 

.40 

2.2443 

2.2079 

2.1053 

1 .9542 

1.7094 

1.2448 

.8597 

.5343 

.2532 

.1361 

.0710 

.0293 

.60 

1.7059 

1 . 6807 

1 . 6094 

1 . 5037 

1.3305 

.9945 

.7073 

.4557 

.2272 

.1263 

.0677 

.0286 

.80 

1.3026 

1 . 2837 

1 . 2302 

1.1510 

1.0212 

.7694 

.5540 

.3645 

.1898 

.1099 

.0613 

.0270 

.90 

1.1338 

1.1166 

1 . 0682 

.9966 

.8800 

.6563 

.4681 

.3060 

.1601 

.0944 

.0542 

.0249 

.95 

1.0554 

1.0387 

.9915 

.9220 

.8093 

.5952 

.4178 

.2680 

.1375 

.0809 

.0470 

.0223 

t  case5 
an\ka 

(2,1) 

.01 

.02 

.05 

.10 

.20 

.50 

1.00 

2.00 

5.00 

10.00 

20.00 

50.00 

.05 

2.6354 

2.5782 

2.4201 

2.1942 

1.8456 

1.2393 

.7907 

.4518 

.1942 

.0988 

.0498 

.0200 

.10 

2.4648 

2.4137 

2.2719 

2.0682 

1.7508 

1.1901 

.7673 

.4427 

.1921 

.0982 

.0496 

.0199 

.20 

2 . 1764 

2.1348 

2.0187 

1 .8503 

1.5838 

1 . 0999 

.7229 

.4248 

.1879 

.0970 

.0492 

.0199 

.40 

1 . 7462 

1.7168 

1.6341 

1.5124 

1.3158 

.9448 

.6411 

.3894 

.1787 

.0941 

.0484 

.0197 

.60 

1 .4359 

1.4135 

1.3504 

1.2570 

1.1045 

.8107 

.5633 

.3516 

.1673 

.0901 

.0471 

.0195 

.80 

1.1921 

1 . 1739 

1.1226 

1 . 0467 

.9225 

.6829 

.4802 

.3053 

.1501 

.0831 

.0446 

.0189 

.90 

1 .0846 

1.0676 

1.0198 

.9491 

.8342 

.6141 

.4300 

.2730 

.1353 

.0760 

.0416 

.0181 

.95 

1 .0331 

1.0164 

.9694 

.9001 

.7878 

.5747 

.3985 

.2504 

.1230 

.0692 

.0383 

.0170 
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A. 4  Vn  interpolation  parameters 


The  sn,  rn,  qn  parameters  from  equation  6.4  are  tabulated  below  with  the  roots  zlt  z2 
of  the  denominator.  Complex  roots  are  again  indicated  by  z\~  z2  —  0. 


# 

Vnnew 

cases  here 

# 

Vnnew 

case=(l ,0) 

# 

an 

sn 

rn 

qn 

zl 

z2 

.05 

15.649078 

6.955146 

2.998574 

.0000 

.0000 

.10 

13.015436 

5.784638 

3.360938 

.0000 

.0000 

.20 

9.049259 

4.021893 

3.705761 

.0000 

.0000 

.21 

8.727861 

3.879049 

3.719977 

.0000 

.0000 

ro 

fO 

8.417930 

3.741302 

3.731308 

.0000 

.0000 

.23 

8.118999 

3.608444 

3.739917 

.0000 

.0000 

.24 

7.830621 

3.480276 

3.745960 

-.4903 

-.5861 

.25 

7.552372 

3.356610 

3.749582 

-.4400 

-.6770 

.26 

7.283846 

3.237265 

3.750921 

-.4159 

-.7428 

.27 

7.024656 

3.122070 

3.750108 

-.3996 

-.8016 

.28 

6.774435 

3.010860 

3.747266 

-.3875 

-.8571 

.29 

6.532829 

2.903479 

3.742512 

-.3781 

-.9109 

.30 

6.299502 

2.799779 

3.735954 

-.3706 

-.9638 

.31 

6.074134 

2.699615 

3.727697 

-.3645 

-1.0164 

.32 

5.856417 

2.602852 

3.717838 

-.3594 

-1.0690 

.33 

5.646058 

2.509359 

3.706472 

-.3552 

-1.1218 

.34 

5.442778 

2.419013 

3.693684 

-.3518 

-1.1752 

.35 

5.246309 

2.331693 

3.679557 

-.3489 

-1.2291 

.36 

5.056393 

2.247286 

3.664171 

-.3466 

-1.2839 

.37 

4.872787 

2.165683 

3.647598 

-.3447 

-1.3396 

.38 

4.695256 

2.086780 

3.629908 

-.3432 

-1.3963 

.39 

4.523574 

2.010477 

3.611168 

-.3421 

-1.4541 

.40 

4.357527 

1.936679 

3.591439 

-.3412 

-1.5132 

.50 

2.965678 

1.318079 

3.351233 

-.3453 

-2.1972 

.60 

1.956834 

.869704 

3.059728 

-.3646 

-3.1535 

.70 

1.219928 

.542190 

2.740393 

-.3959 

-4.6584 

.80 

.679650 

.302067 

2.401634 

-.4408 

-7.5099 

.90 

.284419 

.126408 

2.032365 

-.5081 

-15.5697 

.95 

.129933 

.057748 

1.816074 

-.5606 

-30.8877 

# 

n=  100 

# 

an 

sn 

rn 

qn 

zl 

z2 

.10 

12.931217 

5.747208 

3.368665 

.0000 

.0000 

.50 

2.951357 

1.311714 

3.363324 

-.3433 

-2.2208 

.90 

.284799 

.126578 

2.052668 

-.5028 

-15.7139 

52 


Vnnew 

case=(l , 1) 

an 

sn 

rn 

qn 

zl 

z2 

.05 

7.241101 

7.241101 

2.764161 

.0000 

.0000 

.10 

6.264352 

6.264352 

2.992208 

.0000 

.0000 

.20 

4.700724 

4.700724 

3.258106 

.0000 

.0000 

.30 

3.521932 

3.521932 

3 . 347423 

.0000 

.0000 

.40 

2.616385 

2.616385 

3.320101 

-.4918 

-.7772 

.50 

1.909595 

1.909595 

3.212369 

-.4124 

-1.2698 

.60 

1.350404 

1.350404 

3.045078 

-.3990 

-1.8559 

.70 

.902797 

.902797 

2.827394 

-.4064 

-2.7254 

.80 

. 540847 

.540847 

2.555768 

-.4305 

-4.2950 

.90 

.245451 

.245451 

2.199782 

-.4803 

-8.4819 

.95 

.117948 

.117948 

1.951891 

-.5292 

-16.0195 

Vnnew 

case=  (2,0) 

an 

sn 

rn 

qn 

zl 

z2 

.05 

15.696815 

6.976362 

2.998734 

.0000 

.0000 

.10 

13.170855 

5.853713 

3.369399 

.0000 

.0000 

.20 

9.472064 

4.209806 

3.777178 

.0000 

.0000 

.30 

6.965895 

3.095953 

3.938914 

-.3504 

-.9219 

.40 

5.207810 

2.314582 

3.990995 

-.3042 

-1.4200 

.50 

3.938549 

1.750466 

4.011750 

- . 2846 

-2.0072 

.60 

2.999865 

1.333273 

4.057477 

-.2705 

-2.7727 

.70 

2.291505 

1.018447 

4.189919 

- . 2544 

-3.8596 

.80 

1.748091 

.776929 

4.523079 

-.2302 

-5.5915 

.90 

1.326076 

.589367 

5.404253 

-.1889 

-8.9807 

.95 

1.150844 

.511486 

6.490885 

-.1560 

-12.5343 

Vnnew 

case=  (2,1) 

an 

sn 

rn 

qn 

zl 

z2 

.05 

7.263342 

7.263342 

2.763266 

.0000 

.0000 

.10 

6.340202 

6.340202 

2.993940 

.0000 

.0000 

.20 

4.926668 

4.926668 

3.294619 

.0000 

.0000 

.30 

3.910764 

3.910764 

3.472497 

.0000 

.0000 

.40 

3.156499 

3.156499 

3.595992 

-.4821 

-.6571 

.50 

2.580335 

2.580335 

3.709645 

-.3594 

-1.0782 

.60 

2.128512 

2.128512 

3.850302 

-.3143 

-1.4946 

.70 

1.765214 

1.765214 

4.062959 

-.2802 

-2.0214 

.80 

1.466120 

1.466120 

4.433660 

-.2455 

-2.7786 

.90 

1.214711 

1.214711 

5.221419 

-.2009 

-4.0976 

.95 

1.103268 

1.103268 

6.094277 

-.1693 

-5.3546 
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A. 5  Normalized  crack-mouth  openings  Sn(R)  against  kR 


The  following  data  were  calculated  for  comparison  with  semi-infinite  cracks  from  a 
circular  hole.  Here  the  length  scale  is  R  and  not  a,  hence  spring  stiffness  is  kR. 
The  rows  with  an  —  1.00  were  calculated  for  semi-infinite  cracks  in  a  different  way 
to  finite  cracks  from  the  hole.  The  row  labelled  #dninf  indicates  the  opening  at 
infinity  for  the  semi-infinite  crack  cases. 


#  Normalized  crack  mouth  openings,  dn(l)=  d(x=R)/  R,  d=  2*u,  against  kR  for  semi-infinite 

#  crack  from  a  hole  limit. 

#  holel ,  lbda=  0 


t  an\kr 

.01 

.02 

.05 

.10 

.20 

.50 

1.00 

2.00 

5.00 

10.00 

20.00 

50.00 

.OS 

.0662 

.0661 

.0659 

.0655 

.0647 

.0626 

.0592 

.0535 

.0413 

.0297 

.0189 

.0088 

.10 

.  1273 

.1270 

.1261 

.1246 

.1218 

.1141 

.1031 

.0862 

.0575 

.0366 

.0210 

.0091 

.20 

.2382 

.2371 

.2338 

.2285 

.2185 

.1931 

.1616 

.1216 

.0694 

.0402 

.0218 

.0092 

.40 

.4385 

.4341 

.4214 

.4018 

.3678 

.2935 

.2200 

.1472 

.0747 

.0414 

.0220 

.0092 

.60 

.6558 

.6441 

.6116 

.5646 

.4902 

.3543 

.2458 

.1552 

.0758 

.0416 

.0220 

.0092 

.80 

1.0109 

.9752 

.8833 

.7669 

.6132 

.3964 

.2590 

.1584 

.0761 

.0416 

.0220 

.0092 

.90 

1.4231 

1.3340 

1.1308 

.9155 

.6817 

.4123 

.2628 

.1590 

.0761 

.0415 

.0220 

.0092 

.95 

1.9146 

1 .7166 

1.3373 

1.0121 

.7159 

.4179 
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